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Abstract 

The proximity of the Peccei-Quinn scale to the scale of supersymmetry breaking in models of 
pure gravity mediation hints at a common dynamical origin of these two scales. To demon¬ 
strate how to make such a connection manifest, we embed the Peccei-Quinn mechanism into 
the vector-like model of dynamical supersymmetry breaking a la lYlT. Here, we rely on the 
anomaly-free discrete symmetry required in models of pure gravity mediation to solve the y 
problem to protect the Peccei-Quinn symmetry from the dangerous effect of higher-dimensional 
operators. This results in a rich phenomenology featuring a QCD axion with a decay constant 
of O (l0^°) GeV and mixed WlMP/axion dark matter. In addition, exactly five pairs of extra 
5 and 5* matter multiplets, directly coupled to the supersymmetry breaking sector and with 
masses close to the gravitino mass, m ^/2 ~ 100 TeV, are needed to cancel the anomalies. 
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1 Introduction 

1.1 Common dynamical origin of the PQ and the SUSY breaking scale 

The Peccei-Quinn (PQ) mechanism is widely regarded as the most promising approach to solving 
the strong CP problem in quantum chromodynamics (QCD). Not only does it offer a dynamical 
explanation for the absence of CP violation in strong interactions, it also gives rise to a rich and 
testable phenomenology. In particular, it predicts the existence of a new pseudoscalar particle, 
the axion [^, which could even form dark matter (DM) or at least contribute a sizable fraction to 
it j^. On the other hand, the PQ mechanism on its own is far from constituting a complete theory. 
At the conceptional level, it presumes the existence of an approximate global Abelian symmetry, 
17(1)pq, which is spontaneously broken at some energy scale Apq. In order to account for the tight 
experimental upper bound on the QCD vacuum angle, 6 < 10“^^ [^, this global PQ symmetry 
needs to be of excellent quality, which contradicts with general expectations in the context of 
quantum gravity [^[^. In fact, it is believed that in a consistent quantum theory of gravity all 
global symmetries are necessarily broken by gravitational interactions [^. At low energies, this is 
reflected in a series of higher-dimensional, Planck-suppressed operators which explicitly introduce 
large symmetry-breaking effects. It is, hence, unclear what ensures that the PQ symmetry is 
actually of such high quality and not spoiled by gravitational effects at the Planck scale. Moreover, 
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the origin of the PQ scale Apq remains obscure. In typical axion models 9p0 , the PQ scale usually 
takes a value close to the axion decay constant, Apq ~ O (1 • • • 10) fa, where fa is conventionally 
constrained (based on astrophysical and cosmological observations (^ ) to lie somewhere within 
the so-called axion window, 10® GeV < /a ^ 10^^ GeV. The PQ symmetry is therefore likely to be 
broken at some intermediate scale, which appears surprising, given that one would naively expect 
Apq to be either tied to the electroweak scale Apq ~ 10^ GeV, or, in absence of any other 

scale, to the string or Planck scale, Apq ~ (lO^^ • • • 10^® GeV) / (327r^) ~ 10^^ • • • 10^® GeV 


12 


The first of these two possibilities has, however, been ruled out experimentally a long time ago [P 
while the second one requires fine-tuning of the initial axion misalignment angle to avoid the 
overproduction of axionic DM This leaves one wondering what other new physics might 

possibly come into question as the dynamical origin of the scale Apq. 

The above mentioned problems related to the PQ mechanism appear in a new light as soon 
as supersymmetry (SUSY) is brought into play. Just like the PQ symmetry, SUSY needs to be 
spontaneously broken at some intermediate energy scale, Agusv- This entails the tempting idea to 


suppose that Apq and AgusY are, in fact, determined by the same dynamics 15,^. For instance, 
we may imagine that SUSY is broken dynamically by the interactions in some strongly coupled 
hidden sector. The virtue of such models of dynamical SUSY breaking (DSB) is that, in these 
models, the SUSY breaking scale is generated without any dimensionful input parameters via the 
effect of dimensional transmutation, i.e., it ends up being related to the dynamical scale A of the 
strong interactions, AgusY ~ A. If the same strong dynamics are responsible for the spontaneous 
breaking of the PQ symmetry, we then have 


AgusY ~ A ~ Apq ~ 10^^ • • • 10^^ GeV. 


( 1 ) 


Remarkably enough, this estimate coincides with the range of AgusY values that one typically en¬ 
counters in models of pure gravity mediation (PGM) (l^|19] (for closely related mediation schemes. 


see [M ,21 ). In this framework for the mediation of SUSY breaking to the visible sector, sfermions 
receive large masses of the order of the gravitino mass m^i 2 10^ • • • 10® GeV via the tree-level 
scalar potential in supergravity (SUGRA) [^, while gauginos obtain one loop-suppressed masses 
via anomaly mediation (see also [^). PGM comes with a number of attractive features at 
the cost of a slightly fine-tuned electroweak scale: It easily accounts for a standard model (SM) 
Higgs boson mass of 126 GeV |25| thanks to large stop loop corrections j^; it ensures the unifi¬ 
cation of the SM gauge couplings, so that it may be readily embedded into a grand unified theory 
(GUT); and it is capable of providing a viable candidate for dark matter either in the form of 


winos [19,21,27 , binos 28 or higgsinos 29 . At the same time, it avoids the Polonyi problem 


usually encountered in models of gravity mediation |30|, as SUSY is required to be broken by a 


non-singlet field in PGM; it is free of any cosmological gravitino problem 31 , as the gravitino 


decays way before the onset of big bang nucleosynthesis; and it is in less tension with constraints 


on flavor-changing neutral currents (FGNGs) and CP violation 32 , once again because of the high 
sfermion mass scale. 
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1.2 Unique discrete R symmetry for pure gravity mediation 


Meanwhile, a common idea to protect the PQ symmetry from the dangerous effects of higher¬ 
dimensional operators is to invoke some gauge symmetry which ensures that the PQ symmetry 
survives as an approximate accidental symmetry in the low-energy effective theory [7,33,34. A 
particularly attractive choice in this context, which we have recently examined in more detail 
, is to protect the PQ symmetry by means of a gauged discrete R symmetry, Z^. Such a 


m 


34 


discrete R symmetry may be the remnant of a spontaneously broken continuous symmetry in higher 
dimensions (for instance, in string theory) [^; and it is an often important and sometimes even 
imperative ingredient in SUSY model building: It prevents too rapid proton decay via perilous 


dimension-5 operators 36 ; it forbids a constant term in the superpotential of the order of the 


Planck scale, which would otherwise result in a huge negative cosmological constant |37| ; and 
it may account for the approximate global continuous R symmetry which is required to realize 


stable |38| or meta-stable 39 SUSY-breaking vacua in a large class of DSB models. 


On top of that, a discrete R symmetry automatically suppresses the bilinear Higgs mass in the 


superpotential, W D and hence allows for a solution to the /r problem 40 in the minimal 


supersymmetric standard model (MSSM). In fact, in absence of a bare /r term, a supersymmetric 
mass for the MSSM Higgs fields Hu and needs to be generated in consequence of R symmetry 
breaking. This can, for instance, be done as in the next-to-minimal supersymmetric standard model 
(NMSSM) [^, where one introduces a SM singlet field S that couples to HuH^ and which obtains 
a vacuum expectation value (VEV) {S) = /i during electroweak symmetry breaking (EWSB). A 
more minimal and hence more elegant solution, however, consists in allowing for a Higgs bilinear 
term [42] in the Kahler potential (see also [43| ), 

( 2 ) 


K D cnHuHd , 


where ch is a dimensionless constant of 0(1} and which readily yields a /x term of the order of 
the gravitino mass, /x = 0 x 777 x 3 / 2 , after R symmetry breaking. This solution to the /x problem is 
typically employed in PGM models, where it renders the /x and parameters in the MSSM Higgs 
potential linearly independent, thereby ensuring the successful occurrence of radiative EWSB [44| . 
On the other hand, the term in Eq. Q is only allowed in the Kahler potential as long as the total 
R charge of the Higgs bilinear HuH^ is zero modulo N. In order to assess for which symmetries 
this requirement can be fulfilled, it is important to remember that any discrete symmetry which 
is supposed to be relevant at low energies should be a good gauge symmetry even at the quantum 


level, i.e., it should be free of any gauge anomalies [45 46 . This implies in particular as a min¬ 


imal constraint that the fermions of the MSSM should equally contribute to the weak and color 
anomalies of the discrete R symmetryQM \_Z^-SU(2)l-SU(2)l\ — A [Z^-5C/(3)c'-5'U(3)c'] = 0. 
Eor generation-independent R charges commuting with SU (5) and together with the structure of 


^ Extra matter fields contributing to these two anomalies should appear in complete SU{5) multiplets, so as not 
to disturb the unification of the SM gauge couplings. As such, they, too, should equally contribute to the weak and 


color anomalies of discrete R symmetry, leaving the difference between these two anomalies unaffected 47 . 
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the MSSM superpotential (supplemented by heavy neutrino Majorana mass terms in accord with 


the seesaw mechanism [48] ), this condition then necessitates that (see also 34,49]) 


(N) . 

ruu +rH, = 4 . 


(3) 


Here, we have introduced rx as the general symbol to denote the R charge of the field X and 


(N) 


where = serves as a shorthand notation for equality modulo N. Consequently, only for iV = 4, 
the R charge of the Higgs bilinear vanishes which singles out as the only discrete R symmetry 
consistent with PGM and the generation of the // term via the Higgs Kahler term in Eq. Q. 


1.3 Synopsis: PQ symmetry from dynamical SUSY breaking 

Together, the above observations lead to a remarkably consistent picture in the context of pure 
gravity mediation: While the magnitude of the PQ scale may he determined by the SUSY-breaking 
seale in some appropriate DSB model, the quality of the PQ symmetry may be safeguarded by 
the diserete symmetry required to solve the pL problem. In this paper, we shall demonstrate 
that such a scenario is indeed feasible—we set out to construct an explicit DSB model featuring an 
approximate PQ symmetry that is sufficiently protected by an anomaly-free discrete Zf' symmetry. 
In doing so, we will restrict ourselves to the arguably simplest case, i.e., we will content ourselves 
with presenting a minimal example based on strong SU{2) dynamics breaking SUSY a la lYIT [51] . 
Furthermore, to render the Z^ symmetry anomaly-free, we are lead to introduce a set of new 
SU{5) multipletsj^ which obtain masses via their coupling to the SUSY-breaking sector. Here, the 
phenomenological constraints on our model surprisingly single out a unique number of extra matter 
fields: five pairs of 5 and 5* multiplets with supersymmetric masses close to rn^/ 2 - Such additional 


matter states affect the gaugino mass spectrum in PGM [28,54 and thus play a crucial role in 


determining the composition of the lightest supersymmetric particle (LSP). This has, in turn, 
important consequences for dark matter and SUSY searches at colliders. At the same time, for 
gravitino masses around 100 TeV, our model predicts an axion decay constant fa of O (10^°) GeV 
and an axionic contribution to the total abundance of dark matter of at most O (10 %). Meanwhile, 
the superpartners of the axion, the axino as well as the saxion, receive masses of the order of 


the SUSY breaking scale and hence do not cause any cosmological problems 16 . We therefore 


anticipate our model to offer an appealing solution to the strong CP problem in the context of 
dynamical SUSY breaking, which is consistent with all existing bounds, but which, at the same 
time, can be readily probed in a number of terrestrial experiments and astrophysical observations. 

The rest of the paper is organized as follows: In the next section, we will show how to embed 
the PQ mechanism into the lYIT model of dynamical SUSY breaking. Here, some of the technical 
details regarding the vacuum configuration of the lYIT model after taking into account R symmetry 


^Note that Z 2 is not an actual R symmetry; instead it is equivalent to the non-i? matter parity of the MSSM 


■^We do not consider the possibility of anomaly cancellation via the Green-Schwarz mechanism in string theory 
In such a case, the (and several other discrete R symmetries) could also be rendered anomaly-free solely within 


the MSSM, i.e., without the need for an extra matter sector 53 . 
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breaking have been deferred to Appendix In Sec. we will then comment on the quality of 
the PQ symmetry in dependence of the free parameters of our model and discuss the resulting 
phenomenological constraints. Finally, we will conclude and give a brief outlook in Sec. 


2 Embedding of the PQ mechanism into the lYIT model 

2.1 Field content and superpotential of the SUSY-breaking sector 


The lYIT model represents the simplest vector-like model of dynamical SUSY breaking. In its 
most general formulation, it is based on strongly interacting Sp{N) gauge dynamics, while featuring 
2Nf = 2{N + 1) matter helds T* in the fundamental representation of Sp{N). At energies below 
the dynamical scale A, this theory exhibits a flat quantum moduli space, which is best described 
in terms of the Nf {2Nf — 1) gauge-invariant composite “meson” fields 


Mb'~ , i,j = 1,2,■■■2Nf 

p A ' ' 


(4) 


Here, r/ is a dimensionless coefficient which is supposed to ensure that the meson fields M*-^ are 
canonically normalized at low energies. According to arguments from naive dimensional analysis 
(NBA) [^, it is expected to be of O (Itt); but values as small as ~ vr may perhaps also still be 
admissible. Similar numerical coefficients are expected to appear in a number of places in the low- 
energy description of the lYIT model. For simplicity, we will, however, ignore the possibility that 
these factors could numerically deviate from each other and simply account for the uncertainties 
in all dimensionless couplings of our model by means of a single NBA factor p. 

As shown by Seiberg, the meson fields in Eq. Q are subject to the following quantum mechan¬ 


ically deformed moduli constraint 56 


Pf(Mb) 


(^) 


iV-l-l 


(5) 


where the dynamical scale on the right-hand side of this constraint arises in consequence of non- 
perturbative effects and where Pf (M) denotes the Pfafhan of the antisymmetric meson matrix M, 
[Pf (M)]^ = det (M). A convenient way to implement this constraint when studying the low-energy 
dynamics of the lYIT model is to rewrite it in the form of a dynamical superpotential. 


VEdyn - 


N-l 



( 6 ) 


with the field X acting as a Lagrange multiplier. Unfortunately, the dynamical Kahler potential 
for X is uncalculable and hence the physical interpretation of X is ambiguous. One possibility 
to account for this ambiguity is to define the dimensionless coupling k such that X is always 
canonically normalized, irrespectively of its physical status. More precisely, if a dynamical Kahler 
potential for X should be generated, X would be physical and we would expect k to be some 0(1) 
constant. On the other hand, if X should remain unphysical, we could still introduce X as in 
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Eq. (§, the only difference being that, now, k would formally blow up to inhnity. In what follows, 
we will therefore stick to the dynamical superpotential in Eq. (©, keeping in mind that there are 
two sensible regimes for the coupling k: We should either set k ~ 1 or take the limit k —)> oo. 

Supersymmetry is broken in the lYIT model by lifting the ffat directions in moduli space by 
means of appropriate Yukawa couplings. Let us introduce a singlet field Zij for each and 
couple these singlets to the corresponding “quark” bilinears in the tree-level superpotential, 

= (7) 


At low energies, these Yukawa interactions then give rise to the following effective superpotential. 




1 


A 


eff 

Z Ij 


( 8 ) 


where we have replaced the high-energies Yukawa couplings Ah by their low-energy analogues \ij 
to account for their RGE running from energies high above the dynamical scale to energies below 
the dynamical scale. For all Yukawa couplings \ij being nonzero, the superpotential in Eq. Q 
implies F-term conditions for the meson fields, M = 0, which contradict the deformed moduli 
constraint in Eq. ([^, Pf(M) ^ 0, signaling that SUSY is spontaneously broken. Another way to 
put this is to say that the total effective superpotential at low energies, 


Weff = VEdyn + ~ K 7? Y ( | 


N-1 


Pt{v7 - (f) 


N+l 


+ hij - Zi, , 
2 rj 


(9) 


is of the ORaifeartaigh type and, hence, SUSY is broken via the ORaifeartaigh mechanism 57 


An interesting feature of the lYIT model, which will be of crucial importance to us in the 
following, is that the superpotential in Eq. 0 is invariant under an axial U{1)a symmetry as¬ 
sociated with a global 'k* phase rotation, 'k* —)• exp {iqiO) 'k*, which is anomaly-free under the 


strongly coupled Sp{N) gauge group. In 58 , Domcke et al. have promoted this U{1)a to an 
exact gauge symmetry, U{1)a —> U(1)fi, to point out a possibility how to dynamically generate 
a field-dependent Fayet-Iliopoulos (FI) D-term in field theory. In this paper, we will now identify 
the same global U{1)a as the PQ symmetry, U{1)a —^ f^(l)pQ) and show that it may very well 
serve as a basis for the construction of a viable axion model. 


2.2 SUSY- and PQ symmetry-breaking vacuum at low energies 

In the remainder of this paper, we shall focus on the simplest version of the lYIT model, i.e., SUSY 
breaking via strong Sp{l) = SU{2) dynamics. The extension of our construction to the general 
Sp{N) case is straightforward; it merely requires a bigger notational effort. For Y = 1, we then 
have to deal with Nf = 2 quark flavors, each consisting of a pair of fundamental quarks fields, 
(\k^,\k^) and (T^,^^). Under the U{1)a flavor symmetry, these fields are charged as follows, 

['!''] = ['!'"] = . ['I'l = ['Pi = . (10) 
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where we have chosen the normalization so that the mesons at low energies carry integer charges. 
In fact, relabeling all meson and singlet fields, and according to their U{1)a charges, the 
low-energy effective theory ends up consisting of the following degrees of freedom (DOFs), 

M+ = ^ ^ ^ ^ J^13 ^ J^2 ^ ^14 ^ J^3 ^ J^23 ^ ^4 ^ ^24 ^ 

Z- = Zl 2 , Zj^ = Z 34 , Zq = Zl 3 , Zq = Zl 4 , Zq = Z 23 , Zg = Z 24 . 

In terms of these charge eigenstates, the effective superpotential in Eq. Q now reads 

+ ^ (A+ M+ Z_ + A_ M_ Z+ + Ag Mg Zg) , (12) 

with a = 1,2,3,4, where we have renamed the Yukawa couplings Xij in an obvious manner and 
where the Pfaffian of the meson matrix can now be expanded into the following polynomial, 

Pf (M*^) = M+M_ - Mg^Mg^ -k Mg^Mg^ . (13) 


IFeff ^ KTjX 


Pf(Af«) - ( ^ 


The F-term scalar potential corresponding to this superpotential exhibits a saddle point at the 
origin as well as three local minima, in which the Pfaffian constraint is either approximately satisfied 
by the meson bilinear M+M_, by Mg Mg or by MgMg, with all other meson VEVs vanishing. The 
potential energies of these three vacua respectively scale with the products of the corresponding 
Yukawa couplings, A+A_, AgAg, and AgAg. As we intend to identify the U{1)a flavor symmetry 
with the PQ symmetry, we need to make sure that the U{1)a symmetry is spontaneously broken 
at low energies. This is, however, only achieved once the charged mesons M_|_ and M_ acquire 
nonzero VEVs. In the following, we shall therefore assume that the product A+A_ is (at least 
slightly) smaller than AgAg and AgAg, so that M+M- ~ /rf and Mg = Zg = 0 in the true 
vacuum. Setting all neutral mesons and singlets to zero in Eq. (12), we then obtain the effective 
superpotential describing the fluctuations of M_|_ and M_ around the symmetry-breaking vacuum. 


IVeff K,r]X 


M+M_ - ( - 

T] 


+ - (A+ M+ Z_ + A_ M_ Z+) . 

V 


(14) 


In the limit of rigid SUSY, Z+, Z_ and X turn out to be stabilized at zero (see Appendix [A]) . 
The VEVs of the meson helds M+ and M_ then readily follow from minimizing the F-term scalar 
potential corresponding to the above superpotential. For canonical Kahler potential, we obtain 


M± =e~, X = ^/X^, e = (i_c)V2 
A± r/ 


C= (— 

KT] 


(15) 


where A denotes the positive square root of the geometric mean of A^ and A^ and with e param¬ 
eterizing the suppression of M± w.r.t. the asymptotic expression in the limit k —)• oo (or 0), 


M± = e m 9 , M? = lim M± = 

C^o A± rj 


(16) 


Here, () can be interpreted as a measure for the coupling strength of the Yukawa terms in Eq. (14), 
viz. A, in comparison to the coupling strength of the Lagrange term, viz. kt], which arises in 
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consequence of the deformed moduli constraint. For fixed values of k and r], unitarity and the fact 
that C niust not exceed 1 (so that the meson VEVs do not vanish and the PQ symmetry is, in fact, 
broken) then restrict the Yukawa coupling A to take a value between 0 and Amax — min {dvr, kt]}. 
Consequently, the parameters A, ( and e are allowed to vary within the following ranges, 


0 ^ ^ -^max ) 


0 < c < Cmax = 


A. 


KT] 


(1 Cmax) 


- e 

max; — ‘^mm ^ 


< £ < 1 . 


(17) 


It is important to keep in mind that the above results for the meson VEVs come with a certain 
irreducible uncertainty, given the fact that we are unable to calculate the precise form of the Kahler 


potential below the dynamical scale. In this sense, the exact parameter relations in Eqs. (15), (16) 


and should be taken with a grain of salt, as they may easily receive corrections from the 
noncanonical Kahler potential. On the other hand, we will continue to use the above relations in 
the following analysis for definiteness. After all, they allow for a consistent treatment of our model 
in terms of a well-defined set of parameters with a clear physical interpretation. In other words, 
in the following, we will study the lYIT model for the special case of a canonical Kahler potential, 
guided by the notation that this should represent an important benchmark scenario for the more 
general case —a benchmark scenario that we have well under control. This will, in particular, 
help us to keep track of the various numerical factors in our analysis and provide us with better 
estimates for a handful of prefactors which one would otherwise simply take to be “of 0(1)”. 


Plugging our result for the meson VEVs in Eq. (15) back into the superpotential in Eq. (14), 


we are able to deduce the SUSY-breaking F-terms of the helds and Y, 

A 2 
rf 

The parameter C, can hence also be regarded as a measure for the size of the F-term of the field Y, 


|Fz± I = A (1 - (1 - To , 


\2 a2 

\Fx\ = -2 =C^o 

KT] r]^ 


h? 

Fo = Kr]^. (18) 


-^0 


(19) 


Here, the fact that Fx can be nonzero in the first place indicates (perhaps somewhat surprisingly) 
that, for C > 0, the deformed moduli constraint is actually not exactly satisfied in the true vacuum, 

2 /^n 2 

' ( 20 ) 


M+M_ - ( - 

T] 




where f serves again as a useful measure to parametrize the deviation from the situation in the 
limit K —)> oo (where the moduli constraint is exactly fulfilled in the true vacuum). Moreover, 
given the possibility of nonzero Fx, the field Y may turn out to be the most important one among 


the three SUSY-breaking fields Z_|_, Z_ and Y. According to Eq. (18), we namely have 

1/2 


Fx 


Fz, 


c 


1-C 


( 21 ) 


so that |FA| < |Tz±| for ^ < 1/2 and \Fx\ > |-^z±| for (/ > 1/2 in a nicely symmetric fashion. 
At this point, we should mention that the authors of |58| exclusively focused on the regime of 
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large k and hence small thereby disregarding the possibility that the meson VEVs in Eq. (16) 
could potentially be suppressed. In this paper, we shall, by contrast, refrain from making any such 
assumption regarding the size of ^ and simply keep C as a free parameter in our analysis. This will 
allow us to consistently account for the possibility of suppressed meson VEVs, M± = e , relative 
to the naive expressions which we expect in the limit k —?■ oo. We emphasize that, adapting 
this procedure, we are not only able to capture a possible suppression of the meson VEVs due to a 
nonvanishing F-term for the field V, but—at an effective, phenomenological level—also a possible 
suppression due to the uncalculable strong-coupling effects in the Kahler potential. This is one 
of the main reasons why we decide to stick to our parametrization in terms of C, in the following, 
despite the uncertainties induced by the unknown terms in the dynamical Kahler potential. 


Our results for the singlet F-terms in Eq. (18) now immediately provide us with an expression 
for the SUSY breaking scale Asusy = 


A 2 - ..2 

^'-SUSY 


= A (2 - , 


( 22 ) 


which leads us to yet another interpretation of the parameter From Eq. (22), we infer that 2 — C, 
counts what may be regarded as the effective number of “active” SUSY-breaking fields so 

that ( = 2 — A^gugy More precisely, for ^ = 0, the F-term of the field X vanishes and SUSY is 
solely broken by the F-terms belonging to the fields Z_|_ and Z_. Hence, A'gusY ~ ^ case. On 

the other hand, for C = 1) both Fz_^ and Fz_ are zero. SUSY is then solely broken by the F-term 
of the field X and Aig^gy = 1. Correspondingly, intermediate values of C interpolate between these 
two extrema of A^gugy Furthermore, we are now able to determine the gravitino mass 771 , 3 / 2 , which 
appears as a constant term, Wq = const., in the superpotential upon R symmetry breaking]^ 


W D Wq — raj, 12 Mpj, 


(23) 


with Mpi = (SvrG)”^^^ ~ 2.44 x 10^®GeV denoting the reduced Planck mass. Requiring that the 
true vacuum at low energies be a Minkowski vacuum with (almost) zero cosmological constant, we 
then have to balance the SUSY-breaking contribution to the full SUGRA scalar potential, Agpgy, 
against the constant term, —3/Mpj |lTo|^, induced by R symmetry breaking. This gives 

^ ^ ^iusY ^ A (2 - h? 

^ ^ -s/3Mpi Mpi 

In conclusion, we find that, as anticipated in the introduction, both the SUSY breaking scale as 
well as the gravitino mass turn out to be controlled by the dynamical scale A. In order to attain a 
gravitino mass consistent with PGM, say, 7773/2 100 TeV, we thus need A to be of O (lO^^) GeV. 

In the next section, we shall see how this scale can also be understood as the scale of PQ symmetry 
breaking and, in particular, how it is related to the axion decay constant fa- 

■‘Here, we ignore the VEV of the Kahler potential, Kg, which actnally enters the right-hand side of this relation 
in form of a factor exp (—Ko/Mpi/2). Since Kg Mpj (see Eq. @), this factor is, however, completely negligible. 
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2.3 Identification and decay constant of the axion 


The nonzero VEVs of the charged meson fields M_|_ and M_ in Eq. (15) spontaneously break the 


global U{1)a symmetry of the lYIT superpotential. We identify this flavor symmetry with the 
PQ symmetry, which means that the chiral axion superheld A must correspond to the goldstone 
multiplet of spontaneous U{1)a breaking contained in M_|_ and M_. To make this relation manifest, 
let us expand M_|_ and M_ around their VEVs (which satisfy A+ (M_|_) = A_ (M_)) into a flavor- 
symmetric fluctuation M in the radial direction as well as into a goldstone phase 0, 


M± = ^ 
A± 


A± (M±) + ^M 


.±0 


(25) 


Here, the coupling A/i, denoting the positive square root of the harmonic mean of A^ and A?_, 

- 1/2 


Ah = 


1/1 1 
2 VaI^F 


(26) 


is chosen such that the meson field M, i.e., the “Higgs field” of PQ symmetry breaking, is canoni¬ 
cally normalized. Meanwhile, the phase 0 directly corresponds to the axion multiplet A modulo a 
proper normalization. In order to find the precise relation between 0 and A, we have to examine 
the canonical Kahler potential for M_|_ and M_ in the true vacuum, i.e., for M = 0, 


1 


K = Kocosh (0 + 0t) + A sinh (0 + 0+) = Kq + A (0 + 0t) + -Kq (0 + 0t)^ + 


(27) 


where Kq stands for the VEV of the mesonic Kahler potential and A denotes the difference between 
the two meson VEVs squared]^ 


Ko = ( |M+|2 > + ( |M+|2 >, A = ( |M+|2 > - ( |M+|2 ). 
This leads us to identify the properly normalized axion field A as follows, 


A = kV^Q, Kq = ^{1-C,)(- 


A/i 

" = T = 


1/^ ^ 

2 I A_ A+ 


1 - 1/2 


(28) 


(29) 


Here, the complex scalar (</> -|- i a) /\/2 contained m. A = {(/>, a, a} now consists of the axion a and 
the saxion </>, while the fermionic component of A represents the axino a. Meanwhile, p G [0,1] 


is a convenient measure for the magnitude of the flavor hierarchy in Eq. (14). Eor equal Yukawa 


couplings, A+ = A_, the superpotential in Eq. (14) is invariant under the exchange of and 


and p = 1. On the other hand, as soon as A+ 7 ^ A_, this symmetry is broken and p < ij^ The more 

®Note that A would represent a field-dependent and dynamically generated FI-term for the U(X)a flavor symmetry, 
in case this symmetry was promoted to a gauge symmetry and the meson VEVs were not identical, A+ 7 ^ A- 
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“At high energies and for generic “neutral” Yukawa couplings Ag, the exchange symmetry of the superpotential 


in Eq. (14 1 for A+ = A_ cannot be realized at the level of the fundamental quarks 41*. This renders it an accidental 


rather than an exact symmetry, which is expected to be explicitly broken in the Kahler potential. More generally, 
we stress that, out of the maximal flavor symmetry of the lYIT model, 17(4), we only depend on the PQ symmetry, 
17(4) D S'17(4) D 17(1)pq, to be a reasonably good symmetry. All other global symmetries ought to be explicitly 
broken, in order to avoid massless particles and/or the formation of topological defects in the early universe. 
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we increase the flavor hierarchy in the charged meson sector, the smaller p then becomes—until, 
for A+ ;§> A_ or A+ <C A_, it eventually approaches zero. However, for not too large a hierarchy, 
jlogj^o ("^+/-^)l ^ the parameter p always stays rather close to unity, p > (20/101)^/^ ~ 0.44. 
We will therefore ignore the possibility of a parametrically suppressed value of p in the following 
and simply take it to be some 0(1) constant from now on, i.e., we will work with p ~ 0.3 • • • 1. 


Next, let us determine the axion decay constant fa- As we will discuss in Sec. 3.1, the PQ 


symmetry ends up acquiring a color anomaly due to the presence of additional matter states 
coupling directly to the SUSY-breaking sector. This is good news, since a PQ color anomaly is a 
necessary prerequisite for any implementation of the PQ mechanism; it generates the aGG term 
by means of which the CP-violating 9 term in the effective QCD Lagrangian is eventually canceled. 


.QCD 

h-efr ^%^-Lr 




I . I oig 

^0 




(30) 


Here, Og is the strong coupling constant, G^i, and respectively denote the gluon held strength 
tensor and its dual and Alpg stands for the coefficient of the [/(l)pQ-5t/(3)c-<S't/(3)c’ anomaly. 
The axion decay constant fa is now defined such that these two terms can be combined to yield 


^QCD 


^TV 
faj StT 




(31) 


which tells us that the decay constant fa is basically given by the VEV of the Kahler potential, 

2 e A 


f - ^ ^ 

IaIpqI ° I-ApqI p 


1 (1 _ A)V2 ^ 

V I-ApqI p p 


(32) 


In view of this result, three comments are in order: (i) Irrespectively of the details of how the 
PQ mechanism is actually implemented in a concrete model, naive dimensional analysis leads us 
to expect that fa should be suppressed compared to the PQ breaking scale by one power of the 


NDA parameter p 1^. Our result in Eq. (32) obviously complies with this expectation, (ii) In 


addition to this, the axion decay constant turns out to be further suppressed due to various factors 
beyond the simple NDA estimate, fa ~ A/ 77 . It is also suppressed by the anomaly coefficient 
|AIpq| as well as by the suppression factor e in the meson VEVs. Depending on the size of these 
prefactors, fa may be smaller than the dynamical scale A by one or even more orders of magnitude. 
For e ~ O.I, |AIpq| ~ I • • • 10 and 77 ~ vr • • - dvr, for instance, we typically have fa ~ 10“^ A, so 
that SUSY breaking scales of O (lO^^) GeV result in axion decay constants of O (lO^®) GeV. In 
fact, as we will see in Secs. [3^ and the A and fa values consistent with all phenomenological 
constraints will happen to be of exactly these orders of magnitude, (iii) Similarly to the meson 


VEVs in Eq. (15), the relation between fa and A in Eq. (32) is sensitive to corrections coming from 


the dynamically generated Kahler potential. Imagine, for instance, that the mesons M± obtain a 
quartic Kahler potential due to strong-coupling effects, 

, 2 / , \ 2 


K = M±mI + AK = {1 + c±) M±mI , AK = ±G± Q) , 


(33) 
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where the size of the coefficients C± is unknown and where the c± = Ait'/ parametrize the 

ratio between the noncanonical corrections and the canonical Kahler potential. Large corrections, 
|c±| 0, therefore shift the normalization of the meson helds as well as the VEV of the Kahler 


potential, which in turn modihes the relation between fa and A in Eq. (32). Such large 


corrections do, however, not drastically affect our final conclusions, as all bounds that we are going 


to study in Sec. 3.2 will actually be bounds on A. In our parameter analysis, large corrections 


to the relation in Eq. (32) would therefore only result in different labels along the fa axes in our 


plots; the functional dependences displayed in these plots would still remain the same. 


2.4 Stabilization of the axino and saxion 


Next to the axion a itself, the chiral axion multiplet A in Eq. (29) also contains the superpartners of 


the axion: a two-component Weyl fermion, the axino a, as well as a real scalar, the saxion cj). These 
two particles are potentially produced in large numbers in the early universe by means of various 
thermal and nonthermal processes. This might have significant cosmological implications (^|60] . 
Axino and saxion decays may, for instance, lead to the overproduction of dark matter, inject too 
much entropy into the thermal bath, thereby diluting the primordial baryon asymmetry, or alter 
the predictions of big bang nucleosynthesis. In order to prevent these catastrophic effects from 
taking place, either both the axino and saxion abundances need to be adequately suppressed or 
both species have to decay sufficiently fast. Here, the latter solution is, in particular, realized 
once a and (p are given sufficiently large masses. On the other hand, one can show on rather 
general grounds that, in the supersymmetric limit, the entire axion multiplet A is necessarily 


massless 61 . In the case of unbroken SUSY, the saxion p especially represents a flat direction in 


the scalar potential of the PQ-breaking sector, so that it is prone to cause cosmological problems. 
A successful stabilization of the PQ-breaking fields {M± in our case) can therefore only be achieved 
as long as d and p acquire appropriate soft masses in the course of spontaneous SUSY breaking. 
In the following, we shall show that this is exactly what is happening in our model. 


Eirst of all, let us trade the three SUSY-breaking singlet helds Z_|_, Z_ and X in Eq. (14) for 
the following linear combinations. 


5o = 


^2 = 


1 


(2 - 

1 


(i_^)i/2 ^z+ + z_)-c^/^x 


Si = 


1 


21/2 


iz+ - Z _), 


(34) 


(C/2)^/^ (Z+ + Z_) + 2^/2 (1 - C)^"'^ Y 


(2 - C)^/' 

As we will see shortly, these helds will end up corresponding to the physical mass eigenstates in 


the singlet sector. In terms of these new helds as well as in terms of the helds M and 0 in Eq. (25), 


the effective superpotential of the lYIT model in Eq. (14) can now be rewritten as follows, 

,2 


Wes ~ 


2 2 I 

+ 


m 


2Ai2 

— m she {Fa + M) Si — m 


(2 che Fj - 2 (1 - che) FaM - M^) 


So 


(35) 


r (1 - che) {Fa + M) - i - + 


m 


2r 


.FaM M 


^ 2 , 
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where ch© = cosh© and sh© = sinh0, where we have introduced Fa = Kq = 2 |-^pq| fa as 

an alternative symbol for the normalization of the axion multiplet (see Eq. (29)), where /i = Asusy 


denotes the SUSY breaking scale (see Eq. (22)) and where the parameters m and r are dehned as 


m = Xh 


A 

-, r = 
rj 


c 


2-C 


1/2 


(36) 


Erom the superpotential in Eq. (35), we can calculate the scalar potential for the two scalar 
DOFs contained m A = Fa © D 2“^/^ ((/> + ia), i.e., for the axion a as well as for the saxion (j). As 
discussed in more detail in Appendix [Aj the VEVs of the three singlet fields Sq, Si and S 2 vanish 
in the rigid SUSY limit. Therefore, neglecting any SUGRA effects and setting all singlets to zero, 
the axion scalar potential induced by the spontaneous breaking of SUSY takes the following form. 


V{4>, a) = + m^Fj^ cosh ( V 2 F 




1 




(37) 


Here, we have introduced = 2 m? to denote the saxion mass and used the fact that the four 
parameters //, m. Fa and r are actually not linearly independent; as one may easily check, they 


satisfy the relation (see Eqs. (22), (29) and (36)). The lesson from this scalar 

potential now is twofold. First of all, we note that the scalar potential V{(f>,a) does indeed not 
depend on a, rendering the axion a flat direction. This is, of course, expected, given that the 
field a ought to represent the Nambu-Goldstone boson associated with the spontaneous breaking 
of the U{1)a symmetry by construction. Second, we hnd that the saxion indeed ends up being 
stabilized thanks to the SUSY-breaking dynamics of our model, {(p) = 0. Its mass around the 
origin is controlled by the mass parameter m oc Xh A, which is closely related to the SUSY breaking 
scale and which, moreover, goes to zero as soon as the SUSY-breaking sector decouples from the 
dynamics of PQ symmetry breaking (i.e., as soon as —)• 0 for hxed kt] in Eq. &■ 

In order to study the interactions of the axion multiplet A in the true vacuum, it is therefore 


sufficient to restrict our analysis to the superpotential in Eq. (35) in the limit of small fluctuations 
of the goldstone phase © around zero, 

,2 


lUeff ^ m + - MS2 + 

r 2/i"= 


(^2 _ ^2^ 


(38) 


m 


2 Fa 


1 


2 MAS 1 - -M^S 2 -r {A^ - m2) S 2 


+ 


where we have replaced © by AjFA after expanding in powers of © and where the ellipsis stands for 
operators of dimension 4 and higher. This form of the superpotential provides us with a number of 
useful physical insights: (i) The singlet fields Sq, Si and S 2 indeed parametrize the fluctuations of 
the physical mass eigenstates around the true vacuum. Here, Si turns out to share a Dirac mass 
niA = Tn with the axion held A, while S 2 turns out to share a Dirac mass ruM ^ m/r with the 
“radial” meson held M. Meanwhile, Sq remains massless at tree level, (ii) As is now evident, the 
mass parameter m corresponds to the common Dirac mass of A and Si, while r parametrizes the 
gap in the mass spectrum, i.e., the ratio between the two Dirac masses, r = niA/mM- (hi) Among 
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the three singlet fields, So is the only one with a nonvanishing F-term. We can, thus, identify 
it with the goldstino (or Polonyi) multiplet which is responsible for the spontaneous breaking of 
SUSY via its F-term, \FsJ = /i^. Upon spontaneous SUSY breaking, its fermionic component, the 
goldstino So, is therefore absorbed by the gravitino G (playing the role of its longitudinal DOFs 
thereafter), which is why it eventually acquires a mass rrisQ = Tn^i 2 - At the same time, the scalar 
component of the goldstino multiplet, the sgoldstino sq, is a flat direction of the scalar potential 
at tree level, as it is present in any SUSY-breaking model of the O’Raifeartaigh type. 

At the loop level, the (pseudo)modulus sq C Sq is lifted via radiative corrections. The relevant 


loop diagrams arise from the A'^Sq and M'^Sq Yukawa interactions in Eq. (38) as well as from the 
Yukawa interactions of the X component of the goldstino field, Sq = —rX + • • •, with the neutral 
mesons in the full effective superpotential (i.e., from the Mg Mg A terms in Eq. (|12[)). This results 
in the following contribution to the sgoldstino mass 62 (see Appendix [a| for details) 


mn = 


ujir = 


21 n 2 - 1 
IGvr^ 

1 

21 n 2 - 1 


1+ a;(r)-h ^ a;o(Ao’'^)‘^o(Ao^) 

^1 + ^) 1^(1 + ^')ln(l-r2)- 


m 

— I 771“ 


(39) 


Here, the function lo, which smoothly interpolates between w(0) = 0 and a;(l) = 1, acts as weight 
for the relative importance of the M'^Sq interaction compared to the A'^Sq interaction. In the 
case of a degenerate mass spectrum (i.e., for r = 1), diagrams with virtual M lines in the loop 
yield the same contribution to the sgoldstino mass as diagrams with virtual A lines. On the other 
hand, once the meson field M becomes much heavier than the axion field A (i.e., for r —?■ 0), the 
sgoldstino mass ceases to receive contributions from M loops. In this limit, saxion and axino loops 
then remain as the only source of mass generation for the sgoldstino in the charged meson sector. 
At the same time, the contributions to rrig due to the interaction of Sq with the neutral mesons 
are weighted by wq, evaluated as a function of the Yukawa couplings Ag^ and Ag’^, respectively. 
The full expression for wq is given in Appendix Eor now, we merely remark that, as long as 
Ag > A for all a, also this weight smoothly interpolates between 0 and 1. Here, the maximal value, 
ujQ = 1, is, in particular, attained in the flavor-symmetric limit, i.e., for all Ag being equal to A. 
Furthermore, we note that the contributions to mg induced by the neutral meson loops come with 
prefactors proportional to and A 23 , which are defined as follows. 


Ai4 — 


^■‘ + )(W)'-W)')'' 


-|l/4 


A23 — 


A-+ i ((A§f - (4)7 


-I 1/4 


(40) 


Tn 
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, a similar expression for the loop-induced mass of the pseudoflat direction has been derived (see Eq. (57) 
in this paper). This expression does, however, not feature the weight function lu, as the authors of only work in 
the limit k —> 00 , where cu —^ 0. In this limit, the deformed moduli constraint is fulfilled exactly, the meson field M 
decouples completely and the M^Sq interaction no longer contributes to mo. Similarly, the Mq MqX interactions 


have been neglected in 58 , so that the terms weighted by ojo are missing in this paper. Meanwhile, the calculation 
in 62 does account for the MqMqX interactions. But, as it is also based on the assumption of an exactly fulfilled 
moduli constraint (i.e., on ft —> 00 ), it, too, misses the contribution coming from the M^So coupling. 
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In the limit of equal Yukawa couplings, Aq = A, the prefactors (Ai 4 /A)^ and (A 23 /A)^ therefore also 


reduce to unity, so that, in this limit, the sgoldstino mass in Eq. (39) takes the following form, 

2 In 2 - 1 “ 


AS^ A ^ 


ml = 


IGvr^ 


l+a;(r) + ^^ 


m 

— I m“ 




(41) 


For simplicity and since we do not expect any large flavor hierarchy in the lYIT model, we will 
work with this expression for ml (including the approximation k, r‘^) in the following. 

In addition, working with an even more precise expression for ml (such as, for instance, the one 


in Eq. (39)) would not be of much help for another reason: Unfortunately, next to the perturbative 


Yukawa interactions encoded in the effective superpotential, the sgoldstino mass also receives 
contributions from the effective Kahler potential 


62 


The true sgoldstino mass squared, is 


then given as the sum of ml and some dynamically generated and uncalculable correction. 




= mo + Am^^g . 


(42) 


For large Yukawa couplings, A ~ r^, we expect the uncalculable correction Am'j^ ^ to be of similar 
(but not much greater) importance as the perturbative result m^. We note that this will be the 
more relevant case in the context of our phenomenological study later on (see Sec. [^. For smaller 
Yukawa couplings, A <C ry, on the other hand, we have more confidence in the purely perturbative 
calculation. The upshot of these considerations is that the true sgoldstino mass is, most likely, 
always roughly of the order of the expression in Eq. (41), mg^ ~ mo. On top of that, if we further 
assume Am\^^ to be positive, mo represents a lower bound on the actual sgoldstino mass. 




2 \ 2 
m,o > mo . 


(43) 


This comes in handy, because it allows us to determine a conservative upper bound on the sgold¬ 
stino VEV after taking into account the effect of R symmetry breaking (see Appendix [A]) . Such a 
conservative upper bound on (So) is useful, since it prevents us from underestimating the impact 
of higher dimensional operators on the quality of the PQ symmetry (see Sec. |^. 

(iv) Finally, a few comments on the masses of the remaining bosonic and fermionic DOFs 
contained in Sq, Si, S 2 , M, and A are in order. In the true vacuum and neglecting the effect of 
SUGRA on the VEV of the sgoldstino field Sq, the scalar masses in our model are given as follows. 


~ "To , 


m‘i = 2m? 


mi = 0. 


2 2 
m ± = m , 




m / nx 


2 ^ 
= ,2 


Similarly, we obtain for the fermionic masses in the globally supersymmetric limit 


m-so = 0 , 


2 2 2 
= m 


= 


m 


(44) 


(45) 


The dependence of these different mass eigenvalues on the Yukawa couplings A and k becomes 
more transparent, if we rewrite them as functions of C = A^/ (k^ ry^) (see Eqs. (|15|) and (36)), 


m 




= c, 


m?/r"^ 


= 2-C, 


m ?(1 -|- 

p2i^2^2 


= 2 . 


771' 


’/r^ (1 - r^) 

p2^2^2 


= 2(1-0, (46) 
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These expressions allow us to study the SUSY-breaking and PQ-preserving limit (A —)■ kt]) as 
well as the SUSY-preserving and PQ-breaking limit (A —?■ 0) of our model in a nice fashionj^ For 
X ^ KT], the above masses squared (in units of approach {1,1, 2,0}, while for A —)• 0, they 

turn into {0, 2,2,2}. Here, the massless held in the PQ-preserving limit (the real meson scalar 
m_) is the result of an accidental cancellation in the scalar mass matrix for the special parameter 
choice X = Kt], In the limit A —>■ kt/, the held m_, thus, becomes the second lightest state in the 
lYIT model, the only lighter held being the massless axion a. At the same time, the massless helds 
in the SUSY-preserving limit correspond to the DOFs contained in Sq, Si, and A. We, hence, see 
once again that it is mandatory to break SUSY in order to stabilize the axino as well as the saxion. 

Next, we note that, according to the above results for the mass eigenvalues in our model, the 
physical helds at low energies appear to correspond to four real scalars {4>, a, m+, and rri-), three 
complex scalars (sq, si, and S 2 ), one Weyl fermion (sq) as well as two Dirac fermions ((a, si) and 


(m, 52 ))- In fact, all mass degeneracies in Eqs. (44) and (45) are, however, lifted through SUGRA 
effects—see Appendix where we derive the VEVs of all singlet helds taking into account the 
effect of R symmetry breaking and state the full expressions for all bosonic and fermion masses 


given a nonzero value of (So). The helds in Eqs. (44) and (45) are therefore only quasi-degenerate, 
i.e., they are only degenerate in the rigid SUSY limit. In the full SUGRA case, we have to deal 
instead with ten real scalars, one Weyl fermion and four Majorana fermions. The mass splittings 
among the quasi-complex scalars and quasi-Dirac fermions is then of O ['m^/ 2 ) therefore quite 
large. Last but not least, we mention that, imposing the deformed moduli constraint exactly, i.e., 
in the limit k —)■ 00 , the helds contained in M and S 2 become formally inhnitely heavy. In this 
limit, they are, thus, unphysical and need to be integrated out (see also the discussion in [^). 


The effective superpotential of the lYIT model in Eq. (38) then turns into 


2 

Weff ~ p^So -mASi + ^A'^So . 


2/i^ 

which is nothing but the superpotential studied in (see Eq. (43) therein). 


(47) 


3 Quality of the PQ symmetry and phenomenological constraints 

3.1 Protecting the PQ symmetry by means of an anomaly-free symmetry 

Up to now, we have only discussed the renormalizable interactions among the helds of the lYIT 
model. In the context of SUGRA, we, however, expect gravitational effects at the Planck scale 
to induce further, nonrenormalizable interactions among these helds in the low-energy effective 
theory. A priori, there is no reason why these additional interactions should happen to respect 
the global PQ symmetry that is enjoyed by the lYIT model in the rigid SUSY limit. Instead, 

®As for the sgoldstino, we have ml/{p^ kXAX) = (2 In 2 — 1)(4 — 2 ^ (2 — which turns into 

(2In2 — 1)(2 + p®)/p^A^/(87r^) for A —> Kp (~ 2.3 for A = 47r and p = 1) and into 0 for A —> 0. We also recall that 
lAAS) = {A^/rf) j{ZMpX)!^ (2 — C)- This goes to (A^/p^)/(3M|i)/p^ for A —>■ and to 0 for A —>■ 0. 
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the full effective superpotential as well as the full effective Kahler potential at energies below the 
dynamical scale are expected to contain higher-dimensional operators that explicitly break PQ, 




M3 


Ml, iz: 


M* 


±) 




1 


Mr 




(48) 


with M* denoting an appropriate high-energy cut-off scale close to the Planck scale, M* r\j Mpi. 
These operators result in corrections to the ordinary axion potential in QCD, which causes the 
axion VEV to shift from its desired value, (a) = fa 9, to some displaced value, (a) = fa {9 + A(9), at 
which CP is no longer conserved. In other words, the gravity-induced higher-dimensional operators 
in the effective theory re-introduce a nonzero QCD vacuum angle, A0, (which may easily become 
very large, A0 ;:?> 10“^*^) and, hence, bring us back to the original strong CP problem. 

In order to suppress A9 below the experimental bound, A9 < 10“^°, it is necessary to forbid all 
effective operators that explicitly violate PQ up to some high order. As discussed in the introduc¬ 
tion (see also [^), this is best done by invoking a protective gauge symmetry that eliminates all 
of the relevant dangerous operators from the effective theory. In this paper, we shall, in particular, 
rely on a discrete symmetry, which is well motivated from the perspective of PGM (see our 
discussion in Sec. 1.2). Let us now derive the charge spectrum for such a symmetry in the 
context of the lYIT model and assess which PQ-breaking operators in the effective theory it is 
able for forbid. The characteristics of general Z^ symmetries along with the charge assignment 
for the MSSM fields have already been reviewed in 34 , which is why we will be rather brief in 
what follows. First of all, let us group the fields of the MSSM into complete multiplets of SU{5), 
10 = 5* = (P,£), and 1 = (n'^), since we are only interested in charge assignments 

that are at least compatible with SU (5) unification. The MSSM fields are then charged under the 
discrete Z^ symmetry as follows (for details, see Sec. 2.2.2 and Appendix A in 


34 


(rio, r-5*, ri, m,) = ^ (1, -3, 5, 8, 12) + ^ (1, -3, 5, -2, 2) . 


1 


2 a 


(49) 


5 ' ’ ' ’ ’ ^ 5 

Here, a is an integer that can take any value between 0 and 9. Also, notice that the row vector 


multiplied by 2a/5 on the right-hand side of Eq. (49) encompasses the charges of the MSSM 
fields under the Abelian GUT group U{l)x- This group, sometimes referred to as “fiveness”, 
commutes with SU (5) and can be represented as a linear combination of the weak hypercharge Y 
and the difference between baryon number B and lepton number L, i.e., X = 5 (B—L) — 41". We 
note that the MSSM R charges are not uniquely defined, since the MSSM superpotential (including 
Majorana mass terms for the singlet neutrino fields n^) happens to be invariant under Ziq C U{1)x 
transformations. This ambiguity leaves us with ten different solutions for the MSSM R charges. 

Next, we point out that, solely within the MSSM, the Z^ symmetry turns out to be anomalously 
violated at the quantum level for every possible R charge assignment. This is illustrated by the 
fact that the color as well as the weak anomaly coefficient for the Z^ symmetry are always nonzero, 

= A [zf-SU(3)c-SU(3)c] = 6 + Ng (3 no + rg* - 4) , (50) 

Ai^^ = A [Zf-SU(2)L-SU(2)L] =4 + Ng (3 no + rg* - 4) + (rj^„ + - 2) , 
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where Ng = 3 stands for the number of SM fermion generations. To see that and are 
indeed nonzero (without inserting all possible R charges into Eq. (50) by hand), we have to employ 
Eq. (§ as well as the following relations between the charges rio, r ^*, , and , 

2 rio + rn^ = 2 , r^* + rio + rn^ = 2 , (51) 

which directly result from the MSSM Yukawa interactions. We then hnd for -4.^^ and A^^\ 


A^^^ ^=’ 6-4 Ng -2 . 


( 4 ) 


(52) 


In order to cancel the MSSM contributions to and and, hence, render the symmetry 
anomaly-free, we are therefore led to introduce new matter multiplets that transform under SU{5). 
Here, the easiest possibility is to simply add a certain number of 5 and 5* representations, which 
we shall refer to as Qi and Qi in the following. Given /c G N of such new “quark/antiquark” pairs, 
the total R charge of the extra matter fields needs to satisfy 


= rQ + rQ = 2-Ar, Ar = -f(2 + 4£), £ = 0,1, • • • fc - 1, 


'QQ 


k 


(53) 


so as to make the anomaly coefficients vanish. We, thus, conclude that our axion model predicts 
the existence of new SM-charged fields, without the aid of which we would not able not invoke an 
anomaly-free discrete Z^ symmetry as a protective gauge symmetry for the PQ symmetry. 

Moreover, the new matter helds Qi and Qi come in handy for another reason. In order to allow 
for a successful solution of the strong CP problem, the PQ symmetry in the lYIT sector needs to 
exhibit a color anomaly. This is now easily achieved by coupling Qi and Qi to the SUSY-breaking 
sector. Eor instance, and w.l.o.g., we may assume that all QQ pairs couple to via some 

Planck-suppressed operators. Above and below the dynamical scale, we then respectively have 


2=1 




k 


2=1 


(54) 


where the dimensionless coefficients Cq, and Cq^ are naively expected to be of 0(1) or at most 
as large as dvr. In view of the superpotential terms in Eq. (54), three comments are now in order: 
(i) Since the meson field M_|_ carries PQ charge -|-1, each quark/antiquark pair must carry PQ 
charge —1. Meanwhile, as no MSSM field couples to the SUSY-breaking sector directly, all MSSM 
fields remain uncharged under the PQ symmetry. The total PQ color anomaly, hence, receives 
contributions from the new quark helds only, 

= k {qq + Qq) = -k , qq + QQ =-I ■ (55) 

This renders our axion model a special supersymmetric variant of the KSVZ axion model invented 
by Kim, Shifman, Vainshtein, and Zakharov a long time ago [^. (ii) The superpotential couplings 


in Eq. (54) also act as mass terms for the new quark helds. In fact, upon spontaneous PQ symmetry 
breaking, each quark pair acquires a supersymmetric Dirac mass close to the gravitino mass, 

1/2 


M* 7/^ A+ M* U-cJ 


(56) 
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Therefore, depending on the value of m. 3 / 2 , the extra quark helds may or may not be light enough 
to be within the reach of a future multi-TeV collider experiment. Albeit extremely challenging, 
the discovery of 0{k) new SU{5) multiplets in the vicinity of the gravitino mass would then, of 
course, be a smoking-gun signal of our axion model. 


(iii) Last but not least, we note that the form of the superpotential in Eq. (54) together with 


the lYIT superpotential in Eq. (12) suffices to fix the R charges of all fields in the lYIT sector. 


rM± = ±Ar , rz^ =2± Ar , = 0 , rz^=rx = 2. 


(57) 


Here, we have required that the symmetry be anomaly-free under the strongly coupled SU{2). 
Simply by itself, this implies that the charges of all lYIT quarks must sum to zero: rM+ +?’m_ = 0 
and, thus, rx = 2. Notice that, if this was not the case, the dynamically generated superpotential 
(see Eq. ([^) would explicitly break R symmetry, —)• 0, so that we could no longer rely on 
R symmetry as a tool to constrain the low-energy effective theory. In particular, we would loose 
control over the dynamically generated terms in the superpotential and Kahler potential, which 
might lead to too large a gravitino mass or other unwanted effects. Moreover, we point out that 


our result in Eq. (57) reveals an interesting relation between the charges of the helds in the lYIT 


sector under the local Z^ symmetry and the charges of the same helds under the continuous global 
R symmetry of the lYIT model in the rigid SUSY limit. Under the latter, all meson helds are 
uncharged, while all singlet helds carry charge -|-2. If we denote these global R charges by for 
the individual helds (and if we denote the corresponding PQ charges by q), we arrive at 


r = + q Ar , 


(58) 


which holds for every held in the SUSY-breaking sector. This is to say that, invoking an anomaly- 
free Z^ symmetry in the lYIT sector, we are actually doing nothing else but gauging a discrete 
subgroup of the global x U{1)a_ symmetry of the lYIT model. As we shall demonstrate 

in the next section, this discrete gauge symmetry then allows us to eliminate dangerous higher¬ 
dimensional operators in the effective theory. At the same time, it also hxes the structure of 
the renormalizable interactions in the lYIT model. Without invoking any further symmetry, the 
renormalizable superpotential and Kahler potential could also contain terms such as 

DZl, Zl, ■■■, Z^, Zl, ■■■ . (59) 


So far, we have simply ignored this issue; and now we see that, in general (i.e., for most values of 
the two integers k and i), such terms are automatically forbidden by the Z^ symmetry. 

Einally, we mention that, as a result of the relation in Eq. (58), the SUSY-breaking sector 
on its own turns out not to break the Z^ symmetry—even though the charged meson helds M± 
carry nonzero R charge and obtain large VEVs. Here, the point is that we can always rotate away 
the charges of the charged meson helds by means of a global PQ phase transformation, such that 
r ^ r' = . The strong dynamics of the lYIT sector therefore only break SUSY as well as the 

global PQ symmetry, but leave the gauged Zf- intact, though, Z^ x [/(1)pq —)• Z^. This remnant 
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symmetry is then only broken, —)• Z 2 , by the constant term in the superpotential 

Wq = m 3/2 -^pi) as well as by higher-dimensional operators in the effective theory (such as the 
quark mass term in Eq. (54)). Here, it is interesting to observe that the Z 2 parity that we 
are eventually left with can be identified with the R parity of the MSSM. Our model therefore 
automatically accounts for the origin of R parity in the MSSM. That is, in contrast to many other 
models, it does not require any extension by, say, a gauged B—L symmetry to do so [64]. 


3.2 Constraints on the axion decay constant 

We have not yet uniquely specified all properties of the Z^ symmetry. Our construction still 
exhibits two free parameters: k, the number of extra quark pairs, as well as Ar, the shift in the 
global R charges (see Eqs. (53) and (58)). In the following, we shall now examine for which 


values of these parameters we have a chance of arriving at a viable phenomenology as well as how 
the other parameters of our model (the axion decay constant fa, the gravitino mass m^i 2 , the 
Yukawa coupling A, etc.) are respectively constrained in these different scenarios. 

Eirst of all, we note that, in order to forbid as many PQ-breaking operators as possible, it 


turns out advantageous to choose the integer I in Eq. (53) such that Ar ends up being a fraction 


and not an integer, Ar 0 N. This already rules out scenarios with only one or two extra quark 


pairs from the starl and implies that only the fe-th powers of the fields M± and Z± can appear 


in the effective superpotential as well as the effective Kahler potential. The crucial point here is 
that, for Ar 0 N, the smallest integer multiple of Ar is nothing but k times Ar, 


k Ar = 2 . 


(60) 


For even and odd values of k, the lowest-dimensional PQ-breaking operators in and are 
then respectively given as follows. 


W: 




eff 


D 


m|, z^ 


k even 


K: 




, m ^/2 ^ odd 


eff 


D 


m 3 / 2 M|, 7713/2 ; k even 


m3/2M|, 


k odd 


(61) 


Recall that the mesons M± acquire VEVs of O (A) (see Eq. ([T^), while, in the context of SUGRA, 
the singlets Z+ and Z_ obtain VEVs of O ( 7773 / 2 ) (see Eq. ( |124[ )). Together with 7773/2 ~ A^/Mpi 


(see Eq. (24)), these estimates allow us to assess the order of magnitude of the respectively most 
important corrections to the axion scalar potential, AI 4 , induced by these PQ-breaking operators, 




AV„ 


/ "i3/2\ 

VMpJ 




M, 


PI 5 


■^eff 


AW 


/ "l3/2\ 


/c+2+c 


M, 


PI; 


(62) 


®If the constant term in the superpotential is generated at very high energies (for instance, via gaugino conden¬ 
sation [^), we do not have to worry about any cosmological consequences of R symmetry breaking. In such a case, 
all dangerous topological defects created during R symmetry breaking will simply be inflated away. 

^®For fc = 1, we are unable to forbid tadpole terms for the singlet fields Z± in the renormalizable Kahler potential, 
K D Z±, while for k = 2, we are unable to forbid supersymmetric mass terms for the same fields in the renormalizable 
superpotential, W 3 Zf. These scenarios are, therefore, unfeasible from the very beginning. 
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where c = 1 for even k and c = 0 for odd k. Here, notice that the meson operators require a 


different power counting than the singlet operators (see also Eq. (48)). As the meson fields are, in 
fact, composite fields, ~ (see Eq. (|^), each meson held is actually accompanied by 

one power of the dynamical scale, so that each meson VEV is bound to come with a suppression 
factor of C)(A/Mpi). Therefore, despite the hierarchy between the actual VEVs, {M±) 3> {Z±), 
the effect of the respective meson and singlet operators ends up being comparablej^ 

A A2 


Mpi 


{Mi 


Mpi 


^- 3/2 ~ {Z±) ■ 


(63) 


The main lesson from Eq. (62) is that all PQ-breaking effects induced by the effective Kahler 


potential in Eq. (61) are suppressed compared to the corresponding effects induced by the effective 
superpotential by at least one power of the ratio m 3 / 2 /Mpi. This is perhaps not much of a surprise, 
since the Kahler potential in Eq. (61) is holomorphic in the helds M± and Z±, so that it can 
only contribute to the total scalar potential via pure SUGRA terms. By comparison, the lowest- 
dimensional nonholomorphic terms in are obtained by multiplying the terms in Eq. (61) by 
the R-invariant field products M±M^ and Z±Z^, respectively. These higher-dimensional terms 
then yield corrections to the axion scalar potential which are of the same order of magnitude as 
the corrections induced by the holomorphic terms in K^. In summary, we therefore hnd that the 
PQ-breaking effects stemming from the Kahler potential are always suppressed and that it suffices 
to focus on the PQ-breaking operators contained in the superpotential in the following. 


Let us now be a bit more specific and write down the operators in in Eq. (48) including 


all prefactors, powers of the dynamical scale A, powers of the cut-off scale M*, etc., 




C 


z± 


k! 


1 or 


"l3/2\ 


Z^ 

J 




c 


M± 


1 


{k\f 


rj A 
M* 


M. 


± 


M, 


k-3 


(64) 


with the coefficients Cz^ and Cm± denoting some unknown constants of 0(1) and where the pref¬ 
actor of Z^ is determined by whether the integer k is chosen to be even or odd (see Eq. (48)). The 
most dangerous corrections to the axion scalar potential resulting from this superpotential (deriv¬ 
ing from F-term contributions well as from A-term contributions in SUGRA) are the following. 


AUa = 




(k-l)\ 
Cm± 1 
(kl)^ 


A± 


1 or 

J T] 


k-l 

T 


(65) 


■q K 

Jl 


kurj X*M^ + k\±^Z*^ + (k-2,) m^/2 


M, 


k-l 


k-3 


h.C., 


®®This different power counting in the case of the meson operators represents a distinctive feature of our dynamical 
axion model, which distinguishes it from our earlier axion models presented in [^. In this earlier work, the PQ- 
breaking fields are taken to be elementary fields, M± —> P, P, so that their VEVs do not end up being suppressed 
by a factor of O (A/Mpi). This allows, inter alia, for the possibility of extra quark fields as heavy as the dynamical 
scale, mq^ ~ A, and increases the magnitude of the PQ-breaking terms in the axion potential. In the present paper, 
the mass scale of the new quark fields is, by contrast, tied to the gravitino mass, mg. ~ A^/Mpi ~ " 013 / 2 , and the 
PQ-breaking terms in the axion potential are generally more strongly suppressed. We emphasize that it is these 
differences that explain why we cannot simply use the results for general symmetries obtained in 


34 


and apply 

them to the present scenario in the special case of a symmetry. Instead, a new and dedicated study is necessary. 
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where all chiral helds are understood to represent their scalar components. In order to make the 
dependence of these terms on the axion held value a manifest, we need to expand the charged 


helds M± and Z± around their VEVs (see Eq. (25)). Taking into account the fact that ^ 0 


m 


SUGRA (see Appendix [A]) , we then have for the complex scalars contained in M± and Z±, 


M± = {M±) exp ( ± 


z a 


V2Fa) ’ 


= (^±) exp ± 


i a 


V2Fa) ’ 


( 66 ) 


where we have set all further scalar DOEs contained in M± and Z± to zero. In passing, we also 
mention that, in SUGRA, the scale Fa also receives contributions from the singlet helds Z±, 


Fa = 


Ko = { | M +|2 ) + ( | M _|2 ) + (| Z +|2 ) + (| Z _|2 ). 


(67) 


However, since the VEVs of the singlets are much smaller than the meson VEVs, this is only a 
small correction compared to the globally supersymmetric case. In the following, we shall therefore 


neglect the SUGRA corrections to Fa in Eq. (67) and simply work with the expression in Eq. (29). 


Plugging the expressions in Eq. (66) into the scalar potential in Eq. (65), we hnd that all 


dangerous operators in the axion scalar potential can be brought into the following form, 


AUa A X 


exp ±i —pz - + h.c. 

' V2FaJ 


= V cos 


/ ka \ 

\V2Fa) 


(/c+l—c)/4 


V 




for some appropriate mass scale v that differs from operator to operator. This correction to the 
scalar potential needs to be compared with the instanton-induced scalar axion potential in QCD, 


y(0) 


2 f2 

mala 


1 — COS I 6 — — 
Ja 


, 1/2 


ma = 


'^nf'TT 


l + Z 


600 peV 


lO^" GeV\ 




(69) 


fa V fa 

with rua denoting the axion mass in QGD, which is determined by the mass m^o ~ 135 MeV, the 
vr^ decay constant f^o ~ 92 MeV, the ratio of the up and the down quark mass, z = rriu/md ~ 0.56, 
as well as by the axion decay constant fa 65 . The sum of and AVa is then no longer minimized 


at the CP-conserving held value (a) = fa 6, but rather at (a) = /a (^ + A0), where 


Ae = A0| 


0 sm 


I-^pqI 


e] =Ad< 


0 smf 


A 0 n = 


k 


l-^PQl"ia/a rnlf 


2 

a J a 


(70) 


up to corrections of O {A6q) and where we have used that |Apq| = k (see Eq. (55)). According 
to the experimental bound on the QCD angle, A6q must be smaller than 10~^^, which leads us to 


< 10-2-5 A, ~ 240 keV, A, = {mafaf^ ^ 77 MeV. 


(71) 


The energy scale of the PQ-breaking operators in the axion potential therefore needs to be ex¬ 
tremely suppressed, i.e., it should be even smaller than half the electron mass! Given our estimate 


of the energy scale v in Eq. (68) and taking the gravitino mass to be of 0(100) TeV, this implies 


that scenarios with only A: = 3 or /c = 4 pairs of extra quarks can be safely ruled out, 
(fc+l-c)/4 


/ ^-3/2 \ 

V Mpi ) 


Mpi < 240 keV, 7713/2 ~ 100 TeV, 


k A A^min 


(72) 
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Whether or not k = 5 extra quark pairs are phenomenologically viable is hard to tell in view of 
this rather simplified estimate. The case k = 5, thus, requires a more careful analysis. In fact, as 
we shall see in the following, it turns out that k = 5 new quark pairs are not only viable, but also 
the unique number of new quark pairs that will allow us to satisfy all bounds at the same time. 

In order to constrain scenarios with k > 5 extra quark pairs more precisely, we need to know the 


exact expressions for the energy scale v belonging to the respective terms in AI4 in Eq. (65). These 
expressions simply follow from substituting all fields in Eq. 
well as Eq. (124) in Appendix [A| for our results for {M±) 

2Cz. ^ -x A 1 ^. y ir, yk-i 


with their VEVs (see Eq. (15) as 
(Z±), and (|Ar|), respectively). 


= 


{k-l)\ 


A± 


1 or ^ 

M* ) 


r] 


M, 


k-3 


(73) 


4 2 Cm± 

Vm± = 
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(A:!)^ 


(^kKTj (I A|) (M^z) + kX±^ {Z^) + (/c - 3) mg/a (M±)^ 


(M 






M, 


k-3 


where the additional factors of 2 cancel with the factor 1/2 in Eq. (68). Imposing the requirement 
that these scales be sufficiently suppressed compared to the “axion scale” (see Eqs. (70) and ([7l|)), 

<4 + Vm. < = 10-1° , (74) 

we are then able to derive two /c-dependent upper bounds on the axion decay constant fa, 

■‘'■a \ 


4 I 4 


<A9r^At, 


Z: 
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fa<f?=A^z> iC, P,k,v)F^z 


(k) 




Ml fa<f^M^ =A^;^{C,p,k,p)F'^^> , F^^>=M, 


i(fc) 


(k) 




?(^) 


AjAfp,y'«'=+'» 


(75) 


A /3max 

AC^o —a 


CmM5 


\k) Aik) 

Here, and represent two dimensionless prefactors, the precise values of which depend on 
four crucial parameters of our model: C (i-®-) the Yukawa coupling A in the lYIT superpotential, 
see Eq. (dH)), p (i.e., the flavor hierarchy in the lYIT sector, see Eq. (|29[)), n (i.e., the physical 
status of the Lagrange multiplier field X, see Eq. ([^), and p (i.e., the numerical uncertainty of all 
coupling constants in the effective theory induced by strong-coupling effects, see Eq. @)> 

1 l/(2fc) 


k even: (</, p, k, r/) = 

k odd: (C, p, k, p) = 

Am iC, P, K, p) = 


B 


k-l 


3 (fc-l)/ 2 p, (^p)2A:-3 ^fc-3/2 

(27r)2(''-l) k^Ml p2k (2 - ()5('=-l)/2 

4 X 3^/2 pi ^^^)2(fc-2) ^k-2 _ ^^fc/2-ri 

(27r)2(fc-l) p2fc+3 p2(fc+l) (2 - ()5*^/2-2 

22 A :+1 31/2 p ,2 Qk/2 ^ 1/2 _ ^)fc/ 2 +l 


(76) 
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l/(2k+2) 


l/(2k+2) 


p2{k+l) p2(k-l) pk+2 ^2 _ (^)V2 D 

where we have introduced the symbols B, C, and D for the ease of notation, 

B= (21n2-l)(4-2C + p°) , 

C = l+(l-p")'/' , 

D = 32TTH{2-Cf p^^-^C + p-^ {2-C)C^ +(A:-3) (k77)2 + 


(77) 


24 


























It is illustrative to evaluate the two bounds in Eq. (75) for a few representative parameter 
values. For A: = 5 and k = 6, for instance, and setting Cz = Cm = 1 as well as M* = Mpi, the two 

(k) (k) J - K 

energy scales and ' in Eq. (75) take the following values, 


F^ = Ff^ = ~ 1.1 X 10^^ GeV, ~ 1.3 x 10 ^^ Q^y ^ 


(78) 


which is well above the lower astrophysical bound on the axion decay constant, fa ^ 10® GeV ( 


see 


Sec. 1.1). At the same time, for k = 1, r/ = dvr and assuming identical Yukawa couplings in the 

1 (5, 

l(5) 


lYIT superpotential (i.e., p = 1), and are all of 0(0.1) for almost all values of C, 

0.1. (79) 


~ Af ~ 




Thus, for both scenarios, k = 5 and k = 6, we find that the axion decay constant is typically 
constrained to be at most of O (lO^®) GeV. Given this result, it is worthwhile to to recall that, 
in order to realize gravitino masses of O (100) TeV, we anticipate fa to actually take a value close 
to 10^® GeV (see the discussion below Eq. (32)). Our above estimate of the upper bound on fa is, 
hence, consistent with this expectation—albeit it seems as if fa should be rather close to its upper 
bound in order to allow for the possibility of gravitino masses of O (100) TeV. We will specify 
these statements in the next section, where we will finally present our bounds on fa along with 
the corresponding values of m^i 2 . Before we are able to do so, there is, however, one more issue 
which we need to address. In addition to the upper bounds on fa derived above, the requirement 
of perturbative gauge coupling unification at the GUT scale also results in lower bounds on fa- 
The new quark flavors affect the running of the SM gauge coupling constants. As we take the 
new quark helds to transform in complete SU{5) multiplets, i.e., 5 and 5*, the gauge couplings 
still unify at the same energy scale as in the MSSM, Aqut — 2 x 10^® GeV. Between the new quark 
mass scale, uiq, and the GUT scale, the new quarks, however, contribute to the beta functions of 
the SM gauge coupling constants, which results in a faster running and, hence, a larger value of 
the GUT gauge coupling, acuT = Sgut/ than in the MSSM, ckqut^ — 1/25. Thus, in order 
for our model to be consistent all the way up to the GUT scale, we must require that the SM 
gauge couplings still unify at some perturbative value, i.e., that acuT does not exceed unity. The 
actual value of acuT in our model depends on the details of the MSSM mass spectrum as well 
as on the number and the mass scale of the new quark flavors (i.e., k and ttiq). In particular, it 
increases with k and decreases with mg. For a given MSSM mass spectrum and a fixed value of k, 
the requirement that acuT must remain perturbative can then be translated into a lower bound 
on the mass scale of the new quark flavors 


12 


«GUT = acuT (w-mssm; k, mg) , ogut < 1 => mg > rrfg^ (wmssm; k) . (80) 

^^While in principle this is a correct statement, in practice, we need to be a bit more careful: Because of the 
uncertainty in the low-energy inpnt parameters, the uncertainties in the MSSM mass spectrum, etc., the SM gauge 
couplings do not always unify exactly at Aqut. Instead, the electroweak gauge couplings often unify, ai = 02 , before 
they actually reach the strong gauge coupling as. At the technical level, we therefore have to impose the condition 
that as (and not “aGux”) must remain perturbative, i.e., we require as < 1 at the scale where ai and as unify. 
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Figure 1: Left panel: Lower bound on the masses of the new quark fields, mg’", as function of the gravitino mass 
and the number of extra quark flavors (see Eqs. (82l and (83l). Right panel: New quark mass scale, mg, as well 
as lower bound on the axion decay constant as function of the gravitino mass, /g’^ according to the requirement of 
perturbative gauge coupling unification (i.e., according to the requirement mg > mg'", see Eqs. (87 1 and (89l) for 
k = 5. Here, all other relevant parameters have been chosen as follows: p — I, k = 1, rj = dyr, Cq- = 47r, M* = Mpi. 


In the context of PGM, the MSSM spectrum is basically characterized by two scales: (i) the 
gravitino mass m^j2, which determines the masses of all sfermions as well as of the higgsinos, and 
(ii) the gaugino mass scale ^-3/2; which is related to the gravitino mass via a loop factor in PGM, 
mi/2 ~ ^71,3/2/ (iSvr^), and which determines the masses of the MSSM gauginos. Motivated by the 
perspective of neutralino dark matter, we shall take the gaugino mass scale to be of 0 ( 1 ) TeV and 
treat 7713/2 as a free parameter in the following. Consequently, the lower bound on the new quark 
mass scale, tttq™, then becomes a function of 7773/2 and k only. 

In order to determine tttq™ as a function of these two parameters, we have to solve the renormal¬ 
ization group equations (RGBs) for the SM gauge couplings. We do so numerically and accounting 
for, in total, three different mass thresholds: We set all gaugino masses to 1 TeV, take the masses 
of all other MSSM sparticles to be equal to 7773/2, and assume all new quark flavors to have a 
common mass equal to tttq. Between the Z pole and max {7773/2, tttq}, we simply use the ordinary 
MSSM one-loop beta functions (including the contributions from the new quark pairs), while for 
energies between max {7773/2, 777 q} and Aqut, we perform a two-loop calculation in the DR scheme. 
The idea behind this procedure is that the SM gauge couplings become more sensitive to small 
changes in the beta functions, the larger they are. We should therefore be a bit more careful in 
tracking the running of SM gauge couplings at larger energies than at lower energies, which is why 
we switch from a one-loop analysis to a two-loop analysis beyond the last mass threshold. The 
result of our calculation is shown in the left panel of Fig. which displays 777 q‘“ as a function of 
7773/2 and k. Note that, here, k is treated as a continuous parameter, although, of course, only 
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integer values of k are physically sensible. For k = 5 and k = 6, for instance, and a gravitino mass 
of 100 TeV, we respectively find (setting all other parameters to the same values as in Fig. [^, 


k = 5 


rriQ- ~ 1.8 X 10^ GeV, k = 6 : ~ 2.1 x 10“ GeV. (81) 

The purple line in the left panel of Fig. [^indicates the boundary between two different hierarchy 
schemes that are possible within our set-up. Above the purple line, the new quarks are required 
to be heavier than the MSSM sfermions; below the purple line, they can also be lighter than the 
MSSM sfermions. In passing, we also mention that our numerical result for m-g™ is nicely ht by 
the following analytical expression. 


2.1 X 10“GeV. 


where the powers p and q can be expanded into polynomials in Ak = k — 5.2, 


P = 


5.0-h 2.9 A/c-0.85 AP 
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"i3/2 


(82) 


(83) 


'-Q ^ "^3/2 

Here, k ~ 5.2 corresponds to the k value for which m-g™ ~ ^ 713/2 — 10^ GeV. We emphasize 
that this result for rug™ holds independently of all other details of our axion model. In fact, it 
represents nothing but the universal lower bounds on the masses of k pairs of 5 and 5* multiplets 
imposed by the requirement of perturbative gauge coupling unification for a specific PGM-inspired 
MSSM mass spectrum. For this reason, we believe that it may also be useful in the context of other 
scenarios, where the MSSM particle content is supplemented by further SU{5) representations. 


For given k and rn^i 2 , the constraint on the new quark mass scale in Eq. (82) now implies a 
lower bound on /„. To see this, let us rewrite niQ in Eq. (56) as a function of fa and C- Eqs. (24) 
and (32) allow us to write the gravitino mass as a function of fa and C first, which leads us to 


"^3/2 


4(1-^)V2 VSMpi’ 4 (i_^)1/2M*' 

Requiring mq to be larger than mg“ then provides us with the following lower bound on fa, 

fa > f^Q = ^ ("^3/2, k) M*] , m3/2 = m 3/2 {f^Q^ , c) , 

Cq kp 


Cqk^P^ /, 


(84) 


(85) 


with mg'“ (m 3 / 2 , being given in Eqs. (82) and (83) and with 7713/2 {fa, C) being given in Eq. (84). 

Notice that Eq. ( |85[ ) only represents an implicit dehnition of our lower bound on the axion decay 
ik) 

constant, as /g ^ still appears in the argument of the gravitino mass on the right-hand side. In 
order to evaluate our lower bound on the axion decay constant numerically, it is therefore still 


necessary, for any given set of input parameter values, to solve Eq. (85) self-consistently for /g 


(fc) 


ik) 

Alternatively, we may also trade the C, dependence of /g for a dependence on the gravitino mass. 

ng relation, wh: 

/ 4^/3 7773/2 Mpi 


To do so, we simply have to use the following relation, which immediately follows from Eq. (84), 

-1/2 


C = 1- 


1 + 


I npk^p^f^ 


( 86 ) 
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(k) 

Plugging this relation into Eq. (85) and solving for /g , we obtain the following explicit expression, 


f > fW _ ( 4\/3m3/2 MpA 
^ y ypi K rj j 


1/2 


1 + 


{m3/2,k) 

Cq -v/S 777,3/2 Afpi 



1/2 

1 ) 



(87) 


For k = 5 and k = G, for instance, and taking 7773/2 to be 100 TeV, this bounds evaluates to (again 
setting all other parameters to the same values as in Fig. [^, 


k = 5 : ~ 7.8 X 10^ GeV, k = 6 : ~ 2.1 x 10^^ GeV. 


( 88 ) 


At the same time, Eq. (86) also allows us to rewrite uiq as a function of fa and 777,3/2, 

21 1/4 


CQk^P^ f? 


"IQ = 


M* 


/ 4^/3 7773/2 Mpi 


(89) 


We plot the expressions for fg'^ and mq in Eqs. (87) and (89) in the right panel of Fig. for 
the special case of /c = 5 extra quark pairs. For k = 5 and the values of the gravitino mass that 
we are most interested in, 7773/2 ~ 100 TeV, we again hnd a bound of O (lO^*^) GeV on the axion 
decay constant—which this time is a lower bound and not an upper bound on fa- In summary, 
it therefore seems as if, for k = 5, the axion decay must indeed be of O (lO^*’) GeV, i.e., it must 
neither be much smaller nor much larger than 10^® GeV in order to satisfy all phenomenological 


constraints at the same time (see also our remarks below Eqs. (32) and (79), respectively). In the 
next section, we are now going to specify these statements in a bit more detail. 


3.3 Final results: viable region in parameter space 


The axion decay constant is bounded from above as well as from below (see Eqs. (75) and (85)). 
Thus, in order to asses the viability of our model, we have to search for regions in parameter space 
where not all possible values of fa are ruled out, but which still allow for a viable range for fa, 


/gV/aS/ig. /i;|=min{/f,/£>} . (90) 

Here, the bounds fg'^ and are functions of, in total, nine different parameters, which provides 
us with a lot of freedom when it comes to picking a concrete realization of our axion model. Due 
to this large parametric freedom, the bounds fg^ and can in principle vary over many orders 
of magnitude, so that it turns out impossible to derive a single unique range of values that the 
axion decay constant is confined to. Also, a systematic scan of the nine-dimensional parameter 
scan appears to be difficult (and perhaps also not very revealing). Therefore, we will simply focus 
on certain representative parameter choices in the following, trying to assess what is achievable in 
our model. To do so, let us first recall which nine parameters fg'^ and actually depend on: 

• The number of extra quark pairs, k (see Eq. (|53[)). According to our considerations in Sec. 
the integer k must be A; = 5 or larger (see Eq. ©)• At the same time, increasing the value 
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of k implies an increase in all three bounds on fa- For too many extra quark pairs, the lower 
bound Jq ' will therefore begin to exceed the upper boundary of the phenomenologically 
viable window for the axion decay constant, > 10^^ GeV. For this reason, we will 
restrict ourselves to scenarios with k = 5, k = 6 oi k = 7 extra quark pairs in the following. 


The Yukawa coupling in the lYIT superpotential. A, or alternatively the parameter G [0,1], 
which parametrizes the suppression of the meson VEVs {M±) (see Eqs. (12) and (|15|)). Note 
that C can also always be traded for the gravitino mass m ^/2 via the relation in Eq. ( |86[ ). In 
the following, we will mainly be interested in those values of A (or Q that yield a gravitino 
mass of 100 TeV. This then eliminates the coupling A as a free parameter from our analysis. 


The parameter p G [0,1], which represents a measure of the hierarchy among the Yukawa 


couplings A+ and A_ in the lYIT sector (see Eq. (29)). As evident from Eqs. (75) and (87), 

ik) 

all bounds on fa increase when going to smaller values of p. Here, the lower bound fg 

ik) 

increases, in particular, faster than the upper bound f)^ - In order to maximize the allowed 
region in parameter space, we should therefore choose the parameter p as large as possible, 
p = 1. Interestingly enough, this coincides with the flavor-symmetric limit, A+ = A_, in the 
lYIT sector and, hence, might be regarded as a sensible and well motivated choice. 


• The parameter k, which indicates the physical status of the Lagrange multiplier field X. As 
noted below Eq. K should be either treated as an 0(1) coupling or sent to infinity. In 
the former case (i.e., when the field X is assumed to be physical), larger k values turn out 
to be more advantageous for our purposes, k> 1, because going to larger values of k relaxes 
both the bounds and (the bound is rather insensitive to k). We will therefore 
distinguish between three different cases in the following: k = 1, k = 4, and k —>• oo. 


The NDA parameter p, which captures the numerical uncertainty of all coupling constants 
in the low-energy effective theory induced by strong-coupling effects (see Eq. @)- Larger p 
implies a stronger bound on /„ coming from the meson operators in the superpotential 
(see Eq. (73) and ([76])), which is why we should actually choose p as small as possible, ry ~ tt. 
On the other hand, p is naively expected to be of O (dvr), which is why we will consider two 
cases in the following: p = tt and p = 4 tt. 


The high-energy cut-off scale M* in the PQ-breaking operators in and ( 


'■TO' 


see 


Eq. (48)). As we take these operators to be generated via gravitational interactions, M* 
is expected to be close to the Planck scale. Eor now, we will therefore simply set M* = Mpi. 


The three dimensionless coefficients Cq, Cz±, and Cm± (see Eqs. (54) and (64)), which we 
expect to take values somewhere between 1 and dvr. To maximize the allowed region in 
parameter space, the coefficients Cz± and Cm± should be chosen as small as possible (see 


Eq. (75)), while the coefficient Cq should be chosen as large as possible (see Eq. (85)). We 
will therefore set Czj^ = Cm± = 1 and Cq = dvr in the following. 
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Figure 2: Constraints on the axion decay constant fa and the gravitino mass 7713/2 for the four different parameter 
scenarios specihed in Eqs. ( |91| l, ( |93[ ), and ( |96[ ). The thick solid and dashed green lines represent the upper bounds 
/z ’ and /jj' (see Eq. ([7^), respectively, while the thick purple lines show the lower bound (see Eq. (|85[)). 
The color code and the dashed black lines indicate the suppression of the meson VEVs compared to the asymptotic 
expression (see Eq. ^), i.e., the value of the parameter e = (1 — where = {\/K/rff (see Eq. (EH)- 

Here, large suppression corresponds to a large Yukawa coupling \ ~ kt], while small suppression corresponds to a 
small Yukawa coupling \<^ kt]. The thin solid gray lines indicate the value of AS in integer steps on a logarithmic 


scale, AS = 10 10 


■ . The thick black solid line marks the values of fa and 7773/2 for which AS = 10 
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Let us now study all combinations of the parameters k, k, and r] according to the above list of 
restrictions. Remarkably enough, it turns out that there is actually only one combination, which 
happens to allow for a viable range of fa values for a gravitino mass of 100 TeV! 


Scenario^: k = 5, k = 4, t] = tt => 7713/2 = 100 TeV, fa — 8xlO^GeV. (91) 


In this scenario (referred to as Scenario A in the following), the axion decay constant ends up being 
tightly constrained to a value close to fa — 10^*^ GeV (as expected). In addition, we now see that 
our axion model turns out to yield a unique prediction for the number of extra quark pairs: We have 
to introduce exactly k = 5 pairs of 5 and 5* multiplets—no more, no less. Furthermore, we find 
that the parameter k is required to take a finite value. This means that the Lagrange multiplier 
field X must correspond to a dynamical field. Sending k to infinity (and hence assuming the 
Lagrange multiplier X to be unphysical) is not an option in Scenario A. Finally, if we allow the 
gravitino mass to vary, also our constraints on fa begin to change. This is shown in the upper 
left panel of Fig. which displays the constraints on /„ and 7773/2 in the case of Scenario A. As 
illustrated by this plot, the meson VEVs are always suppressed in Scenario A, M±/M^ ^ 0 - 1 . 
This indicates that, in Scenario A, the Yukawa coupling A is required to be rather large. 

While fa and 7773/2 are found to be tightly constrained in the minimal version of our model (i.e, 
in Scenario A), there are several (almost trivial) possibilities to modify our model, so as to relax 
the bounds on parameter space. For instance, we may assume a different mechanism to generate 


the masses of the new quark pairs than in Eq. (54). So far, we have taken the new quark pairs 
to be coupled to the SUSY-breaking sector via gravitational interactions, i.e., we have taken the 
high-energy cut-off scale in Eq. (54) to be the scale M* ~ Mpi. This, however, does not necessarily 
need to be the case. The new quark pairs might also couple to the SUSY-breaking sector via 
the exchange of GUT messenger fields T' and T' with masses of O (Agut) and transforming as 
fundamentals of both the strongly coupled 517(2) as well as of 517(5), 


D -h M'T'T' 




1 A 


(92) 


In this case, the cut-off scale in the effective quark superpotential is no longer of O (Mpi), but 
rather of O (Agut)- Effectively, such a situation can be accounted for in our analysis by increasing 


the coefficient Cq in Eq. (54) by a factor Mpj/Agut ~ 100. Setting Cq to Cq = 100 x 47r then 
significantly widens the allowed region in parameter space (see the upper right panel of Eig. ([^). 
We shall refer to this scenario as Scenario B, 


Scenario R: k = 5, k = 4, rj = tt , Cq = 100 x dvr. (93) 

Eor 7773/2 = 100 TeV, increasing Cq to such a large value basically removes the lower bound on fa 
coming from the requirement of perturbative gauge coupling unification. The axion decay constant 
then ends up being constrained by the astrophysical bound fa > 10® GeV as well as by 

7773/2 = 100 TeV ^ 10® GeV < /a < 8 x 10® GeV. (94) 
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Here, smaller values of fa require the meson VEVs {M±) to be increasingly suppressed compared 
to the asymptotic expression (see Eq. Too strong a suppression, however, appears 

implausible, both from the standpoint of our explicit calculation as well as according to our general 
expectation regarding the behavior of the strongly coupled lYIT sector at low energies. We there¬ 
fore believe that the axion decay constant has, in general, a tendency of being as large as possible, 
so as to reduce the suppression of the meson VEVs. As for Scenario B, this means that, despite the 
significant relaxation of the lower bound we actually still expect fa to be of O (l0^°) GeV. 

Another trivial possibility to relax the bounds on parameter space is to increase the cut-off 
scale M* by some factor of 0(1 • • • 47r) If we do so starting with Scenario A, it becomes difficult 
to realize gravitino masses of 100 TeV because /g ' increases too drastically. On the other hand, 
combining Scenario B with a larger cut-off scale does provide us with a viable scenario that also 

admits a gravitino mass of 100 TeV. In this case, also k and r] can again be set to different values. 

Scenario C: k = 5 , k = 4, rj = tt , Cq = 100 x (47r)^ , M* = 47r Mpi, (95) 

Scenario D: k = 5 , k = 1, rj = At: , Cq = 100 x (47r)^ , M* = 47 r Mpi . 

Here, we have multiplied Cq by another factor of 47r to keep the ratio Cq/M^ fixed at the same 
value as in Scenario B. The bounds on the fa~'n^ 3/2 parameter space for these two scenarios are 
shown in the two lower panels of Fig. For m ^/2 — 100 TeV, the axion decay constant is again 
bounded by the lower astrophysical bound, fa > 10® GeV, in these two scenarios. At the same 
time, the upper bound on fa now increases by roughly half an order of magnitude, 

Scenario C: /a < 5 x 10^® GeV, Scenario D: /a < 3 x 10^® GeV . (96) 


Guided by the notion that too strong a suppression of the meson VEVs tends to be unrealistic, we 
suppose that also in Scenarios C and D the axion decay constant most likely takes a value close 
to the upper end of the allowed range. That is, once again, we expect fa r\j IQi® GeV. 


Finally, we mention that, relaxing our restrictions on Cq and M* similarly as in Eq. (96), i.e., 
assuming the new quark masses to be generated at the GUT scale and slightly raising the cut-off 
scale M* above the reduced Planck mass Mpy, a number of further interesting scenarios become 
available. For instance, scenarios with k > 5 extra quark pairs now become viable, such as 


Scenario E: k = 6 , k = 1 , rj = Att , Cq = 100 x (47r)^ , M* = dvr Mpi , (97) 

in the case of which the axion decay constant is again required to take a value of O (lO^®) GeV, 
m 3/2 = 100 TeV ^ 7x10® GeV < /„ ~ 5 x 10^® GeV. (98) 

^^Recall that, in our formal calculation (employing a canonical Kahler potential), this is achieved by fine-tuning 
the Yukawa coupling A, so that it increasingly approaches its maximal value Amax = nrj = Att (see Eq. @). 

^■^This may also be desirable from the perspective of flavor-changing neutral currents, which may still be a little bit 
too large for 1713/2 = 100 TeV. Slightly increasing the cut-off scale M* can then help to fully solve the FCNC problem 
also for a gravitino mass of 100TeV, i.e., without the need for going to 7173/2 as large as, say, 1000TeV |32l|66|. 
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Furthermore, we may now also assume that the field X is unphysical and send k to infinity, 

Scenario F: k = 6 , k— )-oo, t] = iw , Cq = 100 x (dvr)^ , M* = dvr Mp\ . (99) 

In this scenario, the axion decay constant is then more or less constrained to a certain value, 

7713/2 = 100 TeV /a ~ 5 X 10^*^ GeV . (100) 

A more systematic study of these (and possibly other) scenarios is left for future work. For now, 
we merely conclude by observing that our model indeed appears to be compatible with all bounds 
in large parts of parameter space. Without any further assumptions, the number of extra pairs is 
fixed to be A; = 5, while the decay constant fa is generally expected to take a value of O (lO^*^) GeV. 


4 Conclusions and Outlook 


In this paper, we have demonstrated how the PQ solution to the strong CP problem might be 
inherently connected to the dynamics of spontaneous SUSY breaking. To give a concrete example 
of our idea, we have embedded the PQ mechanism into the lYIT model of dynamical SUSY 
breaking (i.e. into a strongly coupled SU{2) gauge theory with four matter and six singlet fields), 
which has led us to a particular supersymmetric variant of the KSVZ axion model. As a direct 
consequence of this embedding, we found that the scale of PQ symmetry breaking, Apq, is no 
longer an arbitrary (and somewhat mysterious) input parameter, but rather directly tied to the 
dynamical scale A of the strong interactions in the SUSY-breaking sector, Apq ~ A. As the same 
dynamical scale also determines the scale of SUSY breaking in the lYIT model, AgusY ~ A, a 
PQ scale of O (lO^^ • • • 10^^) GeV then implies a large SUSY breaking scale and, hence, a large 
gravitino mass, 7773/2 ~ 100 TeV. The proposed connection between the dynamics of PQ symmetry 
and SUSY breaking therefore turns out to go very well with the idea of pure gravity mediation. 

Besides that, the notion of pure gravity mediation is also crucial to our axion model for another 
reason: In order to protect the PQ symmetry from the dangerous effect of higher-dimensional oper¬ 
ators induced by gravitational interactions around the Planck scale, one has to invoke a protective 


gauge symmetry—for instance, as proposed in 34 , a discrete R symmetry. Among all possible 
symmetries, pure gravity mediation singles out the special case of a symmetry, which is 
the only discrete R symmetry that allows to generate the MSSM ^ term via a Higgs bilinear term 
in the Kahler potential. As we were able to show, such a discrete Z^ symmetry then manages to 
suppress all PQ-breaking operators in the superpotential and Kahler potential up to a high order, 
thereby ensuring that the PQ symmetry is of sufficiently good quality. Solely within the MSSM, 
however, a discrete Z^ symmetry does not represent a good symmetry, as it is anomalously vio¬ 
lated at the quantum level by SU{2)l and SU{3)c instanton effects. Therefore, in order to render 
the Z^ symmetry anomaly-free, the presence of further SM-charged is requiredFor this reason, 

^®By contrast, solely within the lYIT sector, the Z 4 is anomaly-free and even preserved in the true vacuum. It 
is, therefore, neither broken explicitly nor spontaneously by the strong interactions, which retains R symmetry as a 
useful tool to study the low-energy dynamics of the SUSY-breaking sector (see our discussion related to Eq. Hz!)- 
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we have assumed the existence of k new quark/antiquark pairs (( 5 i,( 5 i) ~ (5,5*) in this paper, 
which obtain masses of the order of the gravitino mass from coupling to the SUSY-breaking sector. 
Provided an appropriate R charge, these new quark fields then cancel the anomalies, which 
puts us in the position to employ the Zf symmetry as a protective gauge symmetry after all. 

The axion model presented in this paper comes with a number of attractive conceptional and 
phenomenological implications. The extra matter fields, for instance, contribute to the running of 
the SM gauge couplings, which increases the value of the GUT gauge coupling constant. Requiring 
that the SM gauge couplings should unify at a perturbative value therefore puts a lower bound 
on the mass scale of the new quark fields. For one thing, this constrains the parameter space of 
our model (i.e., it provides us with a lower bound on the axion decay constant). For another, we 
note that the effect of several new SU{5) multiplets with masses mg ~ ^ 1 , 3/2 on the running of 
the SM gauge couplings might also play the role of a selection criterion in the landscape of string 
vacua. Of course, such an assertion is highly speculative; but we have the feeling that it is worth 
being pointed out nonetheless. Imagine, for instance, that the SM gauge couplings are bound to 
unify at some 0(1} value at the GUT scale. The fact that the new quark fields obtain their masses 
via couplings to the SUSY-breaking sector may then potentially bias the distribution of different 
values of the SUSY breaking scale—maybe the SUSY breaking scale happens to be very large, so 
that m ^/2 ~ 100 TeV, simply because otherwise the SM gauge couplings would run over too long 
a distance between the GUT scale and the new quark mass threshold. This would then alter the 
ratios of the SM gauge couplings at the electroweak scale and, for one reason for another (in the 
context of nuclear and/or atomic physics), maybe exclude the possibility of habitable universes. 

Apart from this perhaps far-fetched speculation, our model also makes a number of predictions 
which are testable in present-day or near-future experiments. First of all, the new quark fields 
may, for instance, be directly detectable in a multi-TeV collider experiment. In this context, it is 
interesting to remark that our model (at least in its simplest form) surprisingly singles out a unique 
number of extra quark pairs; We have to add exactly five pairs of new matter fields. Remarkably 
enough, this leads to a situation where the R charges of the new quark fields, the R charges of the 
fields in the SUSY-breaking sector as well as the R charges of the MSSM fields all look very similar. 
For k = 5 (and only for k = 5), the R charges of all fields in our model turn out to be multiples 
of 1/5. Whether or not this points at something deep remains to be seen; but it is certainly an 
interesting observation. Moreover, we find that a sufficient suppression of all PQ-breaking effects 
typically requires the axion decay constant to take a value not much larger than fa rsj 10^° GeV. 
This has several interesting implications for cosmology. To begin with, let us remark that, in our 
model, the PQ symmetry should be broken before the end of inflation, i.e., the PQ scale needs 
to exceed the inflationary Hubble scale, Apg > Hinf. If this was not the case, dangerous axion 
domain walls (with domain wall number Advv = I-^pqI = ^ > 1) would form during the QCD 
phase transition, dominating the energy density of the universe soon after their production |67| . 
However, if the PQ symmetry is already broken during inflation, we have to pay attention that 
the isocurvature perturbations induced by the axion fluctuation during inflation, 69 ~ Hiiif/(27r), 
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do not violate any of the stringent bounds derived from the precise observations of the cosmic 
microwave background (CMB) 60 . For fa ~ 10^*^ GeV and an initial axion misalignment angle 9 of 
0(1), this constrains the Hubble rate during inflation to a rather small value, ^ 10® GeV. Our 
axion model is therefore only compatible with small-field models of inflation. Or put differently, 
from the perspective of our model, we are led to expect that upcoming CMB experiments will 
unfortunately not be able to see any signs of tensor perturbations in the CMB. That is, if the 
inflationary Hubble rate should indeed be as small as 10® GeV, or even smaller, the tensor-to-ratio 
is at most of 0(10“^®), which is unfortunately out of reach for any planned CMB experiment. 
Besides that, an axion decay constant of 0(10^^) GeV (in combination with 0 ~ 1 and 66 <C 1), 
results in an axionic contribution to the relic density of dark matter of about 0(10%) [^. The 
remaining DM density is then accounted for by weakly interacting massive particles (WIMPs) in 
the form of MSSM neutralinos in our model. For this reason, we are confident that both axion as 
well as WIMP dark matter searches may, in principle, be able to find positive signals. Moreover, 
for fa 10^® GeV, the axion mass lies in the meV range. Such relatively heavy axions could, 
for instance, be searched for in fifth-force experiments searching for axion-mediated long range 


forces |69| or in experiments aiming at measuring the proton electric dipole moment 70 . Given 
the fact that A9 may easily take a value only slightly below the upper bound in our model, 

A6 ~ 10~^^ • • • 10“^^, (see the upper left panel of Fig. such experiments look indeed promising. 

In summary, we therefore conclude that our axion model not only appears to provide an 
interesting link between dynamical SUSY breaking and the PQ mechanism, it also gives rise to 
a rich phenomenology that is going to be tested in current and upcoming experiments. This is 
exciting and hopefully only a first step towards a better understanding of supersymmetry, dark 
matter and the new physics lurking behind the strong CP problem—which, as we believe, should 
certainly star some kind of axion field of dynamical origin as the main protagonist. 
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A Exact vacuum of the lYIT model for canonical Kahler potential 

In this appendix, we compute the exact VEVs of the SUSY-breaking singlet fields in the lYIT 
model, Z-, and A, under the simplifying assumption of a canonical Kahler potential for 
all relevant fields. Of course, this assumption can never hold true exactly, as, in general, strong¬ 
coupling effects will always induce higher-dimensional terms in the effective Kahler potential. Still, 
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we deem a calculation based on a canonical Kahler potential useful for several reasons. First of all, 
we expect it to represent an important benchmark scenario for the more general case—a benchmark 
scenario that we have well under control and that allows us to obtain a better understanding of 
the various parameter dependences in our model. Moreover, since we lack the ability to calculate 
the dynamical corrections to the Kahler potential, there are, in fact, not many alternatives to 
assuming a canonical Kahler potential, if we are interested in more than just some rough order-of- 
magnitude estimates. Starting from a canonical Kahler potential, we are able to derive a consistent 
set of expressions in terms of a number of well-defined parameters. This would, by contrast, not 
be possible, if we also intended to account for the uncertainties related to the effective Kahler 
potential. On top of that, in case the dynamically generated corrections to the Kahler potential 
become smaller and smaller, the expressions that we are going to derive in the following become 
increasingly accurate, approximating the true results in the lYIT model with arbitrary precision. 
All in all, we therefore believe that the simplifying assumption of a canonical Kahler potential— 
while not exactly reflecting the actual situation in the lYIT model—still captures many of the 
aspects that we are interested in and that it is, hence, worth a closer examination. 


The starting point of our analysis is the effective superpotential in Eq. (12) in combination 


with a canonical Kahler potential for all DOFs in the low-energy effective theory. 
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In the true vacuum of the scalar potential corresponding to these input functions (assuming A+A_ 
to be the smallest among the three products A+A_, AqAq, and AqAq), SUSY is broken by the F-term 


of the following linear combination of the fields Z_|_, Z_, and X (see Eqs. (22) and (34)), 
1 


5o = 


(2-C) 


1/2 


{i-C)^^^{z+ + z_)-c^/^x 


|F5j=/i2 = A(2-C)'/'^ 


A2 
^2 


( 102 ) 


Taking into account the spontaneous breaking of R symmetry in the context of SUGRA, this linear 
combination turns out to acquire a nonzero VEV of O ( 1713 / 2 ). To see this, first of all note that R 


symmetry breaking induces a constant term Wq in the superpotential (see Eq. (23)) 


VF Z) Wq — 1773/2 Afpi. 


(103) 


Together with the SUSY-breaking tadpole term for the goldstino field So in the effective super¬ 
potential, VFeff D fM^So (see Eq. (38)), and together with the loop-induced mass for the complex 
sgoldstino sq C Sq in the effective scalar potential, Ves D |so|^ (see Eq. ([4l|)), this constant 
superpotential gives rise to the following total scalar potential for the complex scalar sq, 

Ueff = 177^0 l^ol^ - 21773 / 2 / 7 ^ (so -b So) . (104) 
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We, thus, find that the interplay between the constant superpotential Wq and the tadpole term 
breaks the rotational invariance in the complex sq plane. That is, while the imaginary part of 
So remains stabilized at 0 thanks to the loop-induced mass the real component of sq obtains 
a linear potential proportional to m ^/2 which shifts its VEV from 0 to some value of O 


2 

(Re {so}) = 
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(Im{so}) = 0 
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(105) 
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This result readily translates into expressions for (.^+), {Z-), and {X). All we need to know is the 


inverse of the transformation between the two field bases (Z+, Z_, A) and {Sq, Si, S2) in Eq. (34), 
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Taking into account that the singlets and S 2 do not obtain a nonzero VEV, this leads us to 

= {So) = ^{l-r^f"{So) , W = {So) = -r{So), (107) 

where we have used Eq. (36) to rewrite the ^-dependent coefficients in terms of the parameter r. 


From Eqs. (105) and (107), we now see that all singlet VEVs crucially depend on the loop-induced 
sgoldstino mass m^p. In order to obtain usable expressions for {Z+), {Z-), and (A), we therefore 
need to determine this mass parameter as precisely as possible. This is what we shall do next. 
The sgoldstino mass follows from the one-loop effective Coleman-Weinberg potential, 
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where stands for the total mass matrix of the lYIT model squared, Q denotes an appropriate 
renormalization scale for the low-energy effective theory and c is a constant that is sometimes 
introduced for cosmetic reasons, but which may as well also be simply absorbed into the scale Q. In 
order to evaluate Vcw and determine , we therefore need to compute the entire mass spectrum 
of the lYIT model for a nonzero value of the goldstino field Sq . In the charged meson sector (which 
includes the chiral superfields M±, Z±, and A), the physical mass eigenstates correspond to ten 
real scalars, one Weyl fermion and four Majorana fermions (see also our discussion at the end of 


Sec. 2.4). Here, the bosonic DOFs consist of the axion a, the saxion cj) as well as the real and 
imaginary parts of the complex scalars contained in the goldstino field, Sq , the singlet field Si 
(which shares a Dirac mass with the axion field), the radial meson field, m^, and the singlet 
field S 2 (which shares a Dirac mass with the radial meson field), Meanwhile, the fermionic 
DOFs consist of the goldstino sq as well as four Majorana fermions forming two pairs of quasi-Dirac 
fermions, (a, si) and (m, S 2 ), where d stands for the axino. A straightforward calculation of the 
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bosonic mass matrix in global SUSY and at tree level, accounting for nonzero Sq, then yields 
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while a similar calculation of the fermionic mass matrix provides us with 
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( 110 ) 


At the same time, the neutral meson sector (which includes the chiral superhelds Mq and 
Zq, where a = 1,2, 3,4) also features 5o-dependent mass eigenvalues. The reason for this is the 
coupling of the field X = —r iSq + • • • to the Pfaffian of the complete meson matrix in Eq. (12). In 
fact, the total effective superpotential for the neutral meson fields takes the following form, 

„ 2 \l /2 


lUeff ^ KT] 


rSo-{l- r^) ' 52 - M^M^) + Ag Zg , 


( 111 ) 


which gives rise to eight complex scalars as well as to four Dirac fermions. Here, the masses of the 
four complex scalars, and contained in the neutral meson fields Mq are gives as 
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( 112 ) 


where we have introduced the symbols (T 14 , A 14 , (T 23 , and A 23 , for the ease of notation (see Eq. (40)), 


"14 — 


1 1^/2 ri 

2 ((^ 0 ) +(-^ 0 ) ) > Ai 4 =(A^ + 5 i 4 ) ^ , 5 i 4 = 2 “ (-^o) ) 


1/2 


, (113) 


and similarly for (T 23 , A 23 , and < 123 . By contrast, the masses of the complex scalars Zq contained in 
the neutral singlet fields Zq turn out to be independent of 5o up to corrections of C>( |5o|^). 


- (Aq)^ { -^ + O 


a = 1 , 2 ,3,4. 


(114) 


38 




































Last but not least, the masses of the four Dirac fermions in the nentral meson sector are given as 


= hS)" (^) 
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With the expressions in Eqs. (109), (110), ( 112| ), ( |114 ), and ( |115| ) at our disposal, we now 
know the entire mass spectrum of the lYIT model up to O (| S'o ^). This allows us to evaluate 
the Coleman-Weinberg potential and, hence, determine the sgoldstino mass. Differentiating Vcw 


w.r.t. to So and Sq finally leads us to the following result for rUsQ (see Eq. (108)) 
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where the function / stands for the following combination of logarithms, 

2 1 / 1 \ 2 


fix) = + ^) (1 + 


and where the parameters S 14 , ti 4 , S 23 , and t 23 are defined as follows (see Eq. (113)), 


Ai4 , <^14 Ai4 , ^14 

Sl 4 — - , tl 4 — - , S23 — - , t 23 — - • 

<714 <714 <723 <723 


(118) 


The weight function w in Eq. (116) accounts for the relative importance of loop diagrams with 


internal M lines compared to loop diagrams with internal A lines. Similarly, wq is a measure for 
the relative importance of loop diagrams involving neutral meson fields. It reduces to iv in the 
flavor-symmetric limit, which is characterized by the parameter t going to zero, i.e., <uo(s, 0 ) = uj{s). 
Both weight functions are normalized such that they smoothly interpolate between 0 and 1, 


a;(0) = 0, a;(l) = 1, a;o(s, t = s) = 0 , wo(l)0) = l- 


(119) 


Here, note that t can never exceed s by definition, t < s (see Eqs. (113) and ( 118| )). The two 
functions co and ujq can, moreover, be conveniently approximated by the following polynomials, 


uj{r)^r^, a;o(s,t) s"* - ( - ) r 


( 120 ) 


^®Recall that, for the purposes of this appendix, we are assuming the Kahler potential to be canonical. However, 
in a more realistic context, we would also expect the presence of uncalculable higher-dimensional terms in the Kahler 


potential, which would yield further contributions to the sgoldstino mass (see our discussion related to Eq. (421). 
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which nicely reproduce the exact identities in Eq. (119) as well as the fact that t<;o(s,0) = u;{s). 

Next, let us evaluate the sgoldstino mass in the flavor-symmetric limit, i.e., for all Aq being equal 
to A. This will provide us with a much simpler expression for that approximates the full result 


in Eq. (116) reasonably well as long as there is no large hierarchy among the Yukawa couplings. 


In the flavor-symmetric limit, we are then allowed to perform the following simplifications 

= (523 = 0 , Ai4 = A23 = (714 = (723 = A , S14 = S23 = 1 ^14 = ^23 = 0 . (121) 

The weight function uq therefore simply evaluates twice to unity, so that rug^ turns into. 
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This is our final result for the sgoldstino mass. Plugging it into Eq. (|105[), we find for (Sq) 
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(123) 


(2 In 2 - 1) [1 uj{r) + 4/p6] p® A^ ~ (2 In 2 - 1) (4 - 2 C + p^) A^ ' 

We stress that, allowing for the possibility of a noncanonical Kahler potential, this result receives 
corrections due to dynamically generated contributions to the sgoldstino mass. If we assume these 
contributions to be positive, the above expression for the goldstino VEV can also be understood 
as a conservative upper limit. On the other hand, for negative mass corrections coming from the 
dynamical Kahler potential, the actual value for (So) ends up being larger than the expression in 


Eq. (123). In the main body of this paper, we shall assume that the former of these two possibilities 


is realized. In this case, working with our result in Eq. (123) will then correspond to a conservative 


treatment of the effect of higher-dimensional operators on the quality of the PQ symmetry. 

Einally, we are ready to compute the VEVs of the singlet fields Z+, Z-, and X. Making use 


of the relations in Eq. (107) as well as of our result for (Sq) in Eq. (123), we eventually find 

_^ l /2 (2 _ ^)2 
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(1 ^) i /2 (2 ^)2 16 ^: 
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It is instructive to consider the behavior of these expressions for certain extreme parameter choices. 
For C 0) for instance, all three VEVs become arbitrarily large, {Z±), (|V|) 3> m^/ 2 - The reason 
for this is simply the inverse cubic power of A appearing in all of the above VEVs. The two 


IGtt 


(124) 


prefactors multiplying 167r3/A3 777 , 3/2 in Eq. (124) stay, by contrast, also finite in the limit 0. 
For ^ 1, on the other hand, {Z±) approaches 0, while (|V|) is bounded from below, 
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7773 / 2 . (125) 
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dvr 
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Interestingly enough, we also find that all three VEVs typically turn out to be smaller than 7773/2 
for most values of More precisely, for ktj = dvr, p = 1 and in terms of the coupling A, we obtain 


K 77 = dvr 


(Z±) < 7773/2 for A > 2.0 77, (|V|)< 7773/2 for A >1.577. 


(126) 


only perform these simplifications in order to obtain a more practical (approximate) expression for the 
sgoldstino mass, i.e., we do not assume that any larger global flavor symmetry is actually realized, since this would 
lead, for instance, to problems involving massless particles (see Eq. (1121 as well as the discussion in Footnote j^. 
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This is advantageous from the perspective of our axion model, as it indicates that dangerous 
higher-dimensional operators involving powers of the fields Z_|_, and/or X may not have as 
strong an effect on the quality of the PQ symmetry as one may naively expect (i.e., if one simply 
estimated {Z±) and (|X|) to be some values of 0 ( 7713 / 2 ))- Lastly, we mention that, in the limit 
of the deformed moduli constraint being exactly fulfilled, i.e., for k —>• 00 , the VEV of the singlet 
field X vanishes completely, (|X|) = 0, while {Z±) turns into a simple function of A, 

4 / 5 A /47r\ ^ 

^ {Z^) ^ ^ 0.16 j (-j . (127) 

Here, the fact that X vanishes is consistent with the observation that X becomes inhnitely heavy 
for 7t —)• oo, indicating that, in this limit, it is an unphysical held that needs to be integrated out. 
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